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FOREWORD 


This  report  is  the  eleventh  in  a  series  describing  research  in  re-entry 
phenomenology  performed  by  Phi lco- Ford  Corporation  as  part  of  the  Advanced 
Penetration  Problems  (APP)  project.  The  work  described  herein  was  sponsored 
by  the  U.  S.  Air  Force  and  the  Advanced  Research  Projects  Agency  and  moni¬ 
tored  by  the  Space  And  Missile  Systems  Organization,  USAF,  under  Contract 
No.  F04701-70-C-0130  Modification  P0004. 

Previous  technical  reports  in  this  series  have  presented  the  turbulence 
structure  in  flow  fields  representative  of  those  encountered  in  re-entry 
flight.  Such  experiments  have  included  compressible  axi- symmetric  and  two- 
dimensional  wakes  and  highly  heated  plasma  jets  in  which  the  electron  be¬ 
havior  was  measured  together  with  the  gasdtynamical  turbulence.  New  tech¬ 
niques  and  devices  were  introduced  to  overcome  the  rigors  of  high-speed, 
high- temperature  flows  and  the  greatly  increased  complexities  of  compres¬ 
sible  turbulence.  The  present  report  deals  with  the  mean  flow  from  a 
hypersonic  turbulent  boundary  layer  experiment.  Analog  hot-wire  anemometer 
signals  from  the  sAme  experiment  have  been  stored  on  magnetic  tape,  and  a 
detailed  profile  of  the  turbulent  flow  in  that  boundary  layer  will  be  ana¬ 
lyzed  and  presented  in  a  forthcoming  report. 

In  addition  to  the  support  and  encouragement  of  the  agencies  mentioned 
above,  the  authors  wish  to  acknowledge  the  assistance  of  Lee  Von  Seggern 
and  the  expertise  with  which  Ernest  L.  Doughman  prepared  the  smallest  hot¬ 
wires  (0.00001"  diameter)  ever  to  be  used  in  supersonic  or  hypersonic  flow. 
The  work  further  benefited  greatly  from  previous  work  done'  in  the  same 
wind-tunnel  by  Dr.  John  Laufe*.  of  the  University  of  Southern  California 
and  Dr.  Harry  Ashkenas  of  the  *et  Propulsion  Laboratory  (NASA)  and  their 
coworkers.  '  .  . 


ABSTRACT 


An  experimental  investigation  was  carried  out  to  determine  the  structure  of 
the  turbulent  boundary  layer  on  the  wall  of  the  Jet  Propulsion  Laboratory 
Hypersonic  Hind  Tunnel  at  a  free-stream  Mach  number  of  9.37.  Profiles  of 
flow  properties  were  obtained  from  pitot  pressure,  static  pressure,.  ?nd 
total  temperature  surveys  made  through  the  6  inch  thick  boundary  layer  at 
a  station  160  inches  from  the  nozzle  throat.  Test1,  were  conducted  pri¬ 
marily  for  the  following  conditions:  unit  free-stream  Reynolds  number  of 
127,000  per  inch,  corresponding  to  a  Reynolds  number  based  on  momentum 
thickness  of  36,800,  and  a  wall  to  free-stream  total  temperature  ratio  of 
0.385.  A  cursory  examination  of  the  boundary  layer  was  also  made  at 
unit  stream  Reynolds  number  of  67,000.  The  static  pressure  measureiw  v.ts , 
although  indicating  negligible  streamwise  variation,  revealed  the  existence 
of  a  significant  pressure  gradient  normal  to  the  wall.  With  the  pressure 
at  ^e  wall  approximately  45Z  greater  than  its  free-stream  value.  The  ~ 
data  indicate  that  the  profile  of  total  temperature  ratio  versus  velocity 
ratio  agrees  with  the  linear  Crocco  theory  in  the  sub-layer  region,  but 
in  the  outer  pprtion  of  the  boundary  layer,  closely  follows  the  quadratic 
relation  which  has  been  observed  to  characterize  nozzle  wall  measurements. 
The  data  also  indicate  that  the  sub-layer  is  0.05  to  0.10  inches  thick 
which  corresponds  to  only  several  percent  of  the  total  boundary  layer 
thickness.  Correlation  of  the  velocity  profile  with  the  conventional  in¬ 
compressible  profile  shows  poor  agreement  in  the  sub-layer  and  wall-of-the- 
wake  regions,  with  the  latter  attributed  primarily  to  the  pressure  gradient 
normal  to  the  wall.  However,  the  experimental  skin  friction  coefficient 
agrees  with  the  value  predicted  by  the  Van  Driest  theory  within  several 
percent.  Finally,  the  total  temperature  profile  determined  from  a  hot-wire 
anemometer  traverse  through  the  boundary  layer  is  in  good  agreement  with 
the  thermocouple  measurements. 
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SECTION  I 

s' 

INTRODUCTION 


1.1  PURPOSE 

The  hypersonic  turbulent  boundary  lay°r  has  not  been  studied  r  ru 

l°™r  “Ith  1M  i""’°“a"” *■  r«.tr,  tT‘ 

foM~ikr  r“*ff°r.  *  turbu,e"”  characteristics  of  such  a  layer  which 
sporadic  at  en!r  $  Ta  M*C-h  4'5  Reference  1).  have  received  „&y 

ip  *  1,  ;  w"?  VCt’  a  co"1’let«  ““<*>■  »f.  the  hypersonic  layer, 

specially  of  it,  turbulence,  would  fulfill  several  urveL  needs  I, 

narro-  ”  end  curre"t 

ii »rP“de : \szrti " £riz*i feorrr- As an°th" 

ail  anal'  a'  ^secondary  ohjkctive  has  been  to  record  on  naenetU  tape 

ment tUrsa^hn^  bo1^  “dary  ^  lay  e^may^be  Answer  cd  ^without  Repeat  in^  cht^expciS- 

ment.  Such  a  magnetic- tape  "lihrarv"  ,-u  ■  *  exptn- 

ready  collected  Th  P  llbrar>  of  the  sensor  signals  has  been  al- 

oy  collected.  The  mean  ^average)  signals  necessary  to  describe  the' 

The  «i°WfrP*rti*S  haV'  r‘'d“Ced  a"d  a-  presented  in  thC  report 

ean  t low  properties  are  needed  to  reduce  the  turbulence  data  and  In 

of 3 the ^l"1  ^nt°.COntribu,;ions  the  various  turbulence  modes.  Analysis 

fluctuation  measurements  will  be  presented  in  a  subsequent  report. 

1.2  EXPERIMENT  DESIGN 

The  following  criteria  were  set  in  designing  this  experiment: 

? 

a-  £2£tinuous_fl£w  was  require<^to  allow  tl.e  boundary  conditions 
and  i  temperature)  to  be  well  established  and  *huS  v-H  defined- 

and  to  make  the  turbulence  measurements  possible. 

b'w  boundarY  l?ver.s,  to  increase  uhe  spatial  resolution  of 

the  probes,  to  enable  probing- the  sublayer  and  to  decrease  the  turbulence 
frequencies  to  a  manageable  range.  „  turbulence 

r  uV  ~?-~dirnensional  fj ow ,  to  enable  comparison  with  the  simpler 
turbulent  boundary  layer  theories.  simpler 

.  ,  ^  ‘  ^  h--^h  Reynolds  numbers,  to  enable  the  layer  to  reach  dvnarir 

equilibrium,  presumably  the  asymptotic  fo- r.  *  "  d> na.  ic 

"■V, 

leal  Ine  ~ ^  numbers  higher  than  6,  so  as  to  approach  a  range  of  pract- 
lnterest  and  one  where  turbulence  was  not" previous ly  studied. 
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viva /of  .SSgig^SSSSgf  1°""  than  ab°"r  l?ftr'°f  “  sur* 

lotroducediiI7i,crhelrt5gtrad'>nr  “  *rbltr«?  complexities  wou.d  be 

Regarding  heat  transfer.  the  actical  requirement  is  for  a  cooled  wall 

z’sHe,  t,:sdr:,£t1:lv:r  SIS  “r 

The  above  criteria  could  be  be-t  sari^f.eH  k, 

Continuous  hypersonic  wind  tunnel  at  the  hi.hest'attainable  ^  * 

temperature .  To  obtain  Curbulencef  thicken  the  lav er  and  “d 

a  very  long  flat  plate  would  be  needed  We  cost  of  H  e<|u;llbrate  iC 

would  be  prohibitive  however  rw  C  1  h  3  Plate  ™>del 

sidewall  boundary  tayer  of  »uch  a  tul"te"U“"  W“S  <*r,Vn  ‘°  USi"s  the 
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SECTION  II 

FACILITY  DESCRIPTION 


The  above  requirements  were  well  fulfilled  by  making  use  of  the  two-dimen¬ 
sional  flexible  nozzle  21-inch  Hypersonic  Wind-Tunnel  (HWT)  at  the  Jet 
Propulsion  Laboratory  (JPL)  in  Pasadena,  California.  This  facility  is 
described  in  Reference  2,  and  is  shown  in  general  view  on  Figure  1.  Its 
proximity  \to  Philco-Ford,  flexible  schedule,  and  low  cost  made  this  faci¬ 
lity  very  attractive  for  this  work.  Tne  flow  in  this  tunnel  is  contin¬ 
uous,  the  Mach  number  is  adjustable  from  4.3  to  10.5  and  the  unit  Reynolds 
number  can  approach  150,000  per  inch  (1.8  x  lO^f1).  The  Reynolds  num¬ 
ber  R  *  based  on  layer  thickness  5  were  known  to  be  of  order  600,000  antf\ 
Re?  was  estimated  capable  of  reaching  20,000.  Furthermore,  favorable  ' 
evidence  of  two-dimensionality,  repeatability  and  dynamic  equilibration 
in  the  sidewall  boundary-  layer  of  this  tunnel  had  already  been  collected 
by  Laufer,  Ashkenas  and  Gupta  (Reference  3)  under  conditions  satisfying 
the  criteria  of  Section  I;  this  will  be  discussed  further  in  Section  V. 

It  was  decided  to  make  these  tests  in  the  boundary  layer  growing  over  the 
upper  wall  (ceiling)  of  the  21-inch  high  x  20  inch  wide  tunnel  test  section 
whose  layout  is  seen  in  Figure  2,  and  which  is  a  smooth  internally  water- 
cooled  surface.  The  point  of  measurement  is  at  the  mid-span  position  (10 
inches  from  each  sidewall)  and  lies  160  inches  downstream  of  the  nozzle 
throat.  The  conditions  selected  were: 


Mach  nuober : 

\  = 

9 .4  (nominal ) 

Total  pressure: 

P 

o 

3200  and  1700  cm  Hg 

Total  temperature: 

T 

o 

1000°F 

Reynolds  number: 

Re  = 

127  , (DO  per  inch  -  P 

o 

=  3200 

cm  Hg 

= 

67,000  per  inch  *  P 

o 

=  .1700 

cm  Hg 

Although  most  of  the  tests  were  performed  at  P  =  3200  cm  Hg,  a  single 
boundary  layer  traverse,  with  measurements  made  at  nine  points  ranging 
from  1  to  5  inches  from  the  wall,  was  also  carried  out  at  P  =  1700  cm  Hz. 

During  operation  of  the  Wind  Tunnel,  pressure  and  temperature  data  are 
collected  automatically  using  the  JPL  Wind  Tunnel  Data  Acquisition  System, 
a  computer  controlled  recording  device.  Analog  signals  are  conditioned  ’  5 

at  the  test  site  and  then  transmitted  to  a  central  acquisition  area  where 
they  are  digitized  and  stored  in  the  computer  memory.  The  data  can  be 
sent  back  to  the  site  for  printout  by  a  line  printer  and  display  on  a 
visual  monitoring  channel. 
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FIGURE  1.  VIEW  or  JPI.  hypersonic  wind  tunned  showing  data  acquisition  station 
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SECTION  III 


DIAGNOSTIC  PROBES  AND  PROCEDURES 


The  Local  mean  flow  properties  were  found  by  combining  readings  taken  with 
a  pitot  probe,  a  static  pressure  probe  and  a  total-temperature  probe. 
Basically  these  readings,  when  used  with  aerodynamic  formulas  such  as  the 
ideal  gas  law  and  the  Rayleigh  pitot  relation,  yield  the  local  flow  velo¬ 
city,  density,  longitudinal  velocity,  etc.  Once  these  are^known,  non- 
dimensional  quantities,  integral  properties  and  some  "wall  effects  Cc.r. 
be  computed. 


The  diagnostic  probes  were  mounted  on  an  instrumentation  rake  shown  in 
Figure  3,  with  their  sensing  tips  lying  on  a  straight  line  normal  to  the 
flow  and  parallel  to  the  ceiling.  The  two  outboard  and  the  center  probe 
positions  were  occupied  by  1/8  inch  diameter  pitot  tubes,  with  the  two 
intermediate  positions  occupied  by  a  total-temperature  probe  and  the  hot¬ 
wire  anemometer.  The  probes  were  separated  by  a  distance  of  3  inches. 

As  shown  in  Figure  4,  the  static  pressure  tube  is  attached  under  the  rase 
below  the  center  pitot  probe,  with  its  sensing  orifices  vertically  sepa 
rated  from  the  other  probes  by  5/16  inch. 


To  improve  the  resolution  of  the  pitot  survey  near  the  wall,  several  test* 
were  performed  with  a  miniature  probe  fitted  to  the  tip  of  one  of  the  out¬ 
board  pitot  tubes.  The  miniature  probe  was  comprised  of  a  flattened  tube 
with  an  opening  0.008  inches  high  and  0.060  inches  wide.  In  order  to 
approach  the  wall  without  interference  the  tube  was  bent  S-shaped  so  that 
its  opening  was  vertically  offset  from  the  remaining  probes.  The  hot-wire 
anemometer  was  designed  in  a  similar  fashion.  Figure  4  shows  an  end  view 
of  the  rake,  as  seen  by  the  approaching  flow,  indicating  the  various  on¬ 
sets  which  had  to  be  accounted  for  during  data  reduction. 

The  rake  was  suspended  on  an  actuator  mechanism  by  which  the  distance  \ 
between  the  wall  and  the  probes  could  be  fixed  with  a  precision  of  0.001  . 
During  the  tunnel/run  Y  could  be  changed  continuously  at  arbitrary  speed 
for  a  maximum  stroke  of  8  inches.  Data  were  taken  only  after  the  rake 
h\d  been  at  the  desired  Y  position  for  several  minutes,  so  as  to  allow 
the  probe  readings  to  equilibrate.  During  operation  of  the  tunnel,  a 
difficulty  arose  in  determining  the  Y  =  0  position  of  the  rake.  Eecause 
of  the  pressure  differential  across  the  wall  of  the  wind  tunnel,  the  ceil¬ 
ing  shifted  downward  approximately  0.030  inches.  Although  the  test  section 
was  provided  with  viewing  ports,  the  diameter  of  the  port  was  less  than 
the  height  of  the  test  section.  This  prevented  a  horizontal  line-of-sight 
through  the  test  section  at  Y  positions  near  the  wall,  an J  in  particular, 
at  the  location  of  the  ceiling  itself.  The  following  procedure  was  final¬ 
ly  used  to  define  the  location  of  the  ceiling  during  a  run.  With  air- or . 
the  Y  indication  of  the  rake  with  the  uppermost  pitot  tube  immediately 
adjacent  to  the  ceiling  was  recorded.  The  rake  was  moved  away  from  the 
wall,  tne  tunnel  was  turned  on  and  as  the  rake  was  slowly  moved  upward, 
the  uppermost  pitot  tube  and  the  image  formed  by  its  reflection  from  the 
cei line  were  observed  visually  until  the  two  merged.  This  position  of  the 
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Take ,  which  could  be  determined  within  0.005  inches,  was  defined  as  Y  *  0 
and  all  Y  indications  were  corrected  by  the  difference  between  the  air-on 
and  air-off  zeros. 

v 

J 

The  schematic  of  the  data  recording  system  is  shown  in  Figure  5  (for  clarity, 
the  static -pressure  probe  has  been  omitted).  Pressure,  temperature,  and  Y  « 
position  measurements  are  processed  by  the  Hind  Tunnel  Data  Acquisition  Sys¬ 
tem  and  printed  out  by  the  JPL  line  printer  for  isediate  monitoring  by  the 
test  personnel.  In  addition,  the  center  pitot  pressure,  the  total  temper¬ 
ature  and  the  Y  position  data,  together  with  the  hot-wire  signals,  are  re¬ 
corded  on  magnetic  tape  to  provide  a  permanent  data  file  for  subsequent . 
processing. 

' 

It  is  clear  that  more  diagnostic  probes  were  available  on  the  rake  than 
were  needed  for  determining  the  local  flow  conditions.  The  final  data  re¬ 
duction  utilized  the  data  from  the  center  pitot  tube,*  the  static  pressure 
probe  (corrected  for  Y  offset)  and  the  total  temperature  probe.  Although 
the  latter  is  laterally  offset  from  the  other  two  by  3  inches,  the  evidence 
of  two- dimensionality  presented  iq  Section  V  allowed  the  use  of  this  read¬ 
ing  without  fear  of  error  due  to  lateral  gradient  effects..  The  hot-wire 
probe,  also  offset  by  3  inches  from  the  center  offers  an  independent  mea¬ 
surement  of  local  total  temperature  T0  and  unit  Reynolds  number.  These 
were  used  as  redundant  cross-checks  on  TQ  and  Reynolds  number  (see  Section 
VIII).  The  turbulence  data  collected  with  this  probe  will  be  described 
in  a  later  report. 


i 


*  For  Y  distances  less  than  0.8  inches,  data  from  the  miniature  probe  was  » 
used.  However,  as  shown  later,  the  pressures  recorded  at  the  three  pitot 
probe  positions  were  in  excellent  agreement  for  distances  close  to  the  wall 
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FIGURE  5.  SCHWATIC  OF  DATA  RECORDING  SYSTW 


SECTION  IV 
INSTRUMENTATION 


4.1  PRESSURE  TRANSDUCERS 

Pressure  measurements  were  made  with  Statham  Instruments  absolute  pressure 
transducers  which  utilize  an  unbonded,  fully  active,  strain  gage  bridge  as 
the  sensing  element.  Pitot  or  impact  pressures  were  measured  using  model 
PA208TC-K  transducers,  with  a  range  of  0  -  15  psia,  while  a  model  PA208TC-5 
transducer,  with  a  range  of  0-5  psia,  was  used  to  record  static  pressure. 
Before  each  run  the  transducers  were  subjected  to  full  vacuum  to  adjust 
the  zero  setting  and  then  calibrated  with  a  McLeod  gage  over  the  antici¬ 
pated  pressure  range.  A  check  of  the  zero  setting  was  made  just  prior  to 
and  immediately  after  each  run.  The  pitot  transducers  exhibited  relatively 
stable  characteristics  with  negligible  zero  point  drift.  The  static  pres¬ 
sure  transducer,  however,  was  initially  observed  to  drift.an  amount  cor¬ 
responding  to  approximately  5-107.  of  the  nominal  measured  pressure.  This 
problem  was  alleviated  to  a  large  extent  by  mounting  the  transducer  in  a 
temperature  controlled  water-cooled  housing. 

The  signals  from  the  transducers  were  converted  to  digital  form  for  visual 
display  and  printout,  with  the  pitot  pressure  output  appearing  as  100  counts 
per  cm  Hg  and  the  static  pressure  output  as  200  counts/mm  Hg.  While  the 
error  in  the  pressure  readings  was  approximately  17.,  the  accuracy  of  both 
measurements  was  dictated  by  non-linearity  and  hysteresis  of  the  trans¬ 
ducers  which  was  0.75%'of  full  scale.  This  is  of  the  same  order  as  the 
value  of  the  static  pressure  throughout  the  boundary  layer  and  Of  the  pitot 
pressure  near  the  wall.  The  resolution  of  the  transducer  therefore  may  have 
contributed  to  the  observed  scatter  in  the  static  pressure  data.  Sections 
V  and  VI,  and  to  the  Uncertainties  in  the  velocity  profile.  Sections  IX 
and  X. 

\  , 

Both  measurements  were  examined  for  corrections  due  to  various  probe  effects 
The  static  pressure  measurements  had  to  be  corrected  for  viscous  interaction 
effects  as  described  in  Section  VI.  The  influence  of  thermal  creep  ahd  vi¬ 
brational  degree-of-freedom  effects  on  both  the  static  and  pitot  pressures 
was  examined  and  found  to  be  negligible.  In  addition,  viscous'  »5'  and 
rarefaction  effects'6'  on  the  response  of  the  pitct  pressure  probe  were 
investigated.  l[or  the  large  pitot  probe,  the  corrections  due  to  these  ef- 
-fects  did  not  exceed  17.  and  were  therefore  ignored.  Fpr  the  miniature  pitot 
probe  the  viscous  interaction  effect  was  more ^pronounced,  reading  an  esti-“ 
mated  maximum  of  107.  near  the  wall.  However,  experimental  data  on  viscous 
corrections  is  scarce  and  restricted  to  specific  probe  geometries  which 
•differ  considerably  from  the  probe  used  in  the  present  tests.  For  this 
reason,  and  since  the  possible  error  did  not  exceed  107.,  no  correction  was 
applied  and  the  miniature  probe  data  was  used  as  recorded. 


4.2  TOTAL-TEMPERATURE  1R0BE 

**  ^  ^ 

As  explained  previously,  the  total  temperature  probe  consisted  of  a  type  J 
iron- cons tantan  thermocouple.  The  thermocouple  leads  were  inserted  through 
and  cemented  in  a  3/32  inch  diameter  alumina  tube  which  was  then  installed 
in  the  1/8  inch  O.D.  copper  tube  located  on  the  instrumentation  rake.  The 
alumina  body  was  notched  and  the  adjacent  tube  wall  grooved  so  that  the 
thermocouple  could  be  held  in  position  by  a  spring  slip  inserted  over  the 
tube,^  The  thermocouple  junction  was  positioned  7/16  inches  from  the  open¬ 
ing  of  the  copper  tube  and  three  .031  inch  diameter  vent  holes  were  drilled 
through  the  tube  wall  at  this  same  location. 

i  ^ 

The  signal  from  the  thermocouple  was  also  converted  to  digital  form  for 
display  and  printout  with  80  counts  corresponding  to  1  mV.  r The  accut^cy  ‘ 
of  the  probe  measurement,  however,  was  dictated  by  the  precision  with  which 
the  tunnel  supply  temperature  could  be  controlled.  During  l  run,  the  sup¬ 
ply  temperature  varied  typically  by  +  2  to  3°C  and  Similar  variations  were 
reflected  in  the  measured  total  temperature. 


The  total  temperature  probe  was  calibrated  in  the  free- stream  of  the  HWT 
at  Mach  9.5  and  6  and  800°K  and  at  Machs4  and  stagnation  temperatures  rang- 
lng  from  300  to  810  K.  The  high  Mach  nunAer  data  was  supplied  by  the  JPL/ 
USC  investigators'^)  while  the  Ma ch  4  data  was  obtained  during  the  present 
study  because  of  concern  that  the  recovery  factor  of  the  prob?  might  de¬ 
pend  on  Mach  number  and  temperature.  Calibration  results  are  shown  “in 
Figure  6  where  the  probe  recovery  factor,  defined  as  the  ratio  of  measured 
temperature  to  actual  stagnation  teaperature  is  plotted  against  Reynolds' 
number  based  on  probe  diameter  and  stagnation  temperature.  It  is  seen 
that  while  the  Mach  6  and  Mach  9.5  data  Coincide,  the  Mach  4  data  shows  a 
lower  recovery  factor  for  a  given  Reynolds  number.  The  Mach  4  data,  how¬ 
ever,  shows  no  significant  separation  due  to  stagnation  temperature.  In 
order  to  delineate  regions  within  the. boundary  layer  where  the  calibration 
data  would  apply,  the  Y  locations,  together  with  the  associated  Mach  num¬ 
ber  and  stagnation  temperature,  corresponding  to  specific  Reynolds  numbers 
have  been  indicated  in  Figure  6.  It  is  clear  Chat  the  high  Mach  number 
data  is  applicable  to  the  outer  half  of  the  boundary  layer.  In  addition 
the  possibie  error  introduced  by  the  separation  in  the  data  at  low  Reynolds 
numbers  is  at  most  several  percent.  Furthermore,  the  Mach  4  data  indicates 
a  possible  tendency  to  merge  with  the  higt  Mach  nbmber  data  at  Reynolds 
numbers  between  100  and  1000.  Consequently,  it  was  decided  to  represent 
the  probe  calibrations  by  the  single  curve,  determined  from  a  least  squares 
Hf-  fo  .he  Mach  6  and  9.5  data  shown  in  Figure  6. 


4.3  HOT-WIRE  ANEMOMETER 


0  *  .  * 

4.3.1  PROBE  CONSTRUCTION  4  ,,  I 

The  anemometer  probe,  shown  in  Figure  7,  consisted  of  a  0.00001  inch  dia¬ 
meter  Pt  -  10Z  Rh  wire  mounted  across  the  tips  ofJ  tWo  sharp-pointed  prongs. 
The  wires  were  mdunted  with  sufficient  slack  to  eliminate  spurious  signals 
due  to  strain-gage"  effects.  For  the  wires  used  in  these  experiments  the 
aspect  ratio  ranged  from  200  to  400. 
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.  i  nf  a  1/32  inch  diameter  x  li  inches  long  alumina 

The  probe  body  consisted  of  3/  alumina  cylinder,  used  to 

cylinder.  A  1/32  inch  diameter  x  k.  the  larger  cylinder 

support  the,  hot-wire  prongs  ng-M***  °oth  cylinders  were  provided 
at  an  angle  of  approximately  rhr®  »h  which  the  hot-wire  leads  were 

with  two  holes  along  their  e“8  wire  encased  in  high  temperature 

inserted.  The  leads  were  30  g  g  PP  because  of  its  low  electrical 

fiberglass  ‘■“‘“i°°;i«^ueriog,  and  bending'c.pabi  lity  without  b.Mk- 

:"ut^:;--r^i,^e!°us.d  A  the 

soldered  on  one  end  of  the  copper  leads.  T"'  *  *  Pyr'werao  cement, 

tubes  were  held  together  with  torning  Class  ^u^ln5uUtor>  and  cures 

easU^wTt^an^orygeo-acetylene  torch,  leaving  a  smooth  glassy  finish. 

.  —  .  „  -./->?  inch  diameter  alumina  tube  was  first  cut 

To  assentole  the  probe,  the  3/3  shaped  grooves  were  cut  in  both 

to  length  with  a  diamond  blade  saw  ™d  *  *h  Sing  «ca .  Next,  the  1/32 
ends  to  provide  the  pyroceram  cement8  real ter^b  ^8  After  the  nickel 

inch  diameter  alumina  tube  was  cu  wire  it  was  inserted  through 

support  wire  was  silver  soldered  to  ^  cured  with  a  torch,  with 

the  alumina  tubes.  2“"th”  tip^of  the  nickel  prongs  free  and 

cleam  ‘of'ce^t*  ^support  wires  were  then  prepared  and  adjusted  for 
mounting  the  hot  wire.  e 

„  Wwr  wire  consists  of  bonding  the  wire  on  the 
The  technique  for  mounting  the  ho  resinate  and  diluted  wi^th 

►  „iri,  ortirf  which  had  been  dissolved  in  a  resinate  t-' 

supports  with  gold,  wh  .  oven  at  1000°F  leaving  a  thin 

toluene.  *  The  resinate  is  ba  e  y  0nly  one  appiication  of  the 

gold  film  approximately  D  micro- inches  -  diameter  wires,  although 

gold  solution  is  required  lor  ~unt  ng  £or  a’good  bond, 

for  larger  wires,  as  many  as  six  co  g  2-mil  Wona5t0n  wire 

In  the  process  of  mounting  the  ware,  a  s““,Un®5"  “  ,,  held  in  a  pair 

with  a  core  of  0.00001  ^  off>  exposing 

“““JJ:  ,ARi"^e  T^.“robe  is  held  in  a  vise  in  the  object  plane 
the  it  -  107.  Rn  wir  v  order  to  observe  the  mounting  cf 

°f  a  60  pOW" xttar2rgageCw“eris  used  to  apply  the  gold  solution  to  the 
the  wire#  Next  a  zu  6a6c  T*  uv  rarpfullv  niani* 

'tips  and  along  the  outside  l'“Sth  of  t  *  P™®*'  tha  ^uCside  of  one  prong, 
pointing  the  tweezers,  the  hot  wire  is  s  outside  of  the 

across  its  tip  and  over  to  the  tip  of  and  up  along  “  u  solutlo„.  gy 

second  prong.  Hie  hot-wire  ^  ”ej|  wI11  break  where  it  joins  the 

wiggling  the  tweezers,  the  ft  107. “ s“"'„tad>  the  probe  is  placed  in  an 

";“*:t°ab“oiteioo?!CFe,^  i!rL  to  three  minutes  the  solution  is  baked  away 

leaving  a  gold-bonded  hot-wire. 

4.3.2  PROBE  CALIBRATION 

Following  the  probe  assembly,  the 

wire  was  —red  in  wire  resistance,  R, 

Five  temperatures  in  this  range  w  due  to  these  currents  is 

-irs-s.'r.ar^r^  rsi  —  -  —  - 
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zero  current  resistance  was  obtained  by  extrapolation  yielding  a  curve  of 
R  versus  temperature,  T.  The  zero  temperature  resistance  Rr  and  the  first 

coefficient  of  resistivity,  Vr,  were  then  found  from  the  resistivity  formula 

\ 

R  *  R  .1  +  y  (T  -  T  )„  (1) 

v  r  y  r 

o  *  -3  -3 

where  Tf  *  0  C>  Typical  values  of  "tx  ranged  from  1.1  x  10  to  1.8  x  10 

ohm/°C,  with  the  majority  corresponding  to  1.4  -  1.5  x  10“  .„ 

i 

It  was  originally  intended  to  .-ins  tall  each  wire  in  the  Philco-Ford  Mach  3 
tunnel  to  check  for  "strain-gage”  effects  and  to  carry  out  p  preliminary 
flow  calibration  prior  to  its  use  in  the  JPL/HWT.  However,  since  the  flow 
environment  in  the  Philco-Ford  tunnel  is  more  severe  than  the  HWT  (larger 
dynamic M>ressure),  wire  survivability  was  very  poor.  Consequently,  this 
phase  was  omitted  and  flow  calibration  of  each  wire  was  carried  out  in  the 
HWT  immediately  prior  to  performing  turbulence  measurements.  The  objective 
of  this  step  was  to  obtain  the  variation  of  the  Kusselt  number  at  zero 
overheat,  Nuot  and  of  the  wire  recovery  factor,  Tj,  with  Reynolds  number, 

Re0.  The  subscript  o  refers  to  stagnation  conditions,  that  is. 


(2) 


(3) 

(4) 


where  q  is  the  heat  transfer  from  the  wire,  d  is  the  wire  diameter,  Tw  is 

the  wire  temperature,  T a*,  is  the  adiabatic  wire  temperature,  k^,  is  the 

thermal  conductivity  of  air,  s  aud  u  are  the  local  free-stream  density  and 

velocity  of  air  respectively,  aud  ^  is  the  viscosity  of  air  evaluated  at  T  . 

o 

The  calibration  was  performed  by  locating  the  wire  in  the  wind  tunnel  free- 
stream  where  the  flow  conditions  are  well  known  and  varying  the  Reynolds 
number  by  setting  the  tunnel  total  pressure  at  specific  values  between  700 
and  3200  cm  Hg.  At  each  pressure  the  derivative  dR/ci^  was  found  from  the 
slope  of  the  R  versus  i^  curve  obtained  by  operating  the  wire  at  8  overheat 
currents.  The  Nusselt  number  at  zero  overheat,  NuQ,  was  then  given  by 


*  From  initial  calibrations  up  to  500°C  it  was  found  that  the  second  co¬ 
efficient  of  resistivity  was  zero.  It  was  possible  therefore  to  restrict 
subsequent  calibrations  to  200°C,  thereby  facilitating  the  calibration 
process . 
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4.4  WIND  TUNNEL  SUPPLY  CONDITIOfB 

The  wind  tunnel  supply  pressure,  which  was  assumed  equal  to  the  r,cf  c 

■ode  1  ^ PA73 1TC3 1M- 3 so"  FT”'*’  W3S  “asu«d  using  a  Statham  Instruments  ^ 
■odel  PA731TC- 1M-350  strain  gage  transducer,  vi  tl.^-ramie  of  0  innn  < 

a  sensitivity  of  0.011  „f  ,„Ie  and  ,  ’uZU^/^uli  of'fiT .  !*’ 

c™„tePPlLtrrr:tU”  “S  “aS“”d  “Ith  a  shielded^chromel-aluiael1  thermo- " 

P  •  Similar  to  the  other  Measurements  the  siena  1*  fmm  pu-  * 

““If— 

Sr  “ — —  5-**=  svsa 


HEAT  TRANSFER  CHARACTERISTICS  OF  TYPICAL  HOT  WIRE  PROBES 


FIGURE  9.  RECOVERY  TEMPERATURE  CHARACTERISTICS  OF  TYPICAL 
NOT  WIRE  PROBES 


SECTION  V 


* 


PRELIMINARY  TEST  RESULTS 


The  data  described  in  this  section  have  been  included  to  demonstrate  the 
nature  of  the  boundary  layer  flow,  the  two-dimensionality  of  the  boundary 
layer,. and  the  reproducibility  of  the  measurements.  In  addition,  the  mean 
flow  properties  calculated  from  these  preliminary  measurements  were  essen¬ 
tial  for  the  proper  design  of  the  hot-wire  anemometer.  Data  is  presented 
which  was  acquired  over  a  period  of  six  months.  Most  of  the  measurements 
were  collected  by  the  USC/JPL  team 'of  investigators. 

Figure  10  shows  the  pitot  profiles  obtained  from  the  two  outboard  pitot 
probes  at  stations  112  and  160  inches  from  the  nozzle  throat.  The  data 
shown  was  reproduced  directly  from  pressure  traces  obtained  with  a  dual 
channel  x-y  plotter.  Because  of  the  finite  size  of  the  recording  pens  the 
Z  *  -6  pressure  trace  was  physically  offset  by  the  amount  indicated  at  the' 
left-hand- side  of  the  Figure  10  ( Ly  =  0.10  inches).  Therefore  the  corrected 
trace  for  this  transducer  is  obtained  by  translating  all  data  points  to  the 
left  by  the  amount  Ly .  In  spite  of  this  shift  in  the  Z  *  -6  curve,  the 
spanwise  agreement  in  the  pressure  profiles  at  both  stations  is  quite  good, 
indicating  that  the  flow  can  indeed  be  considered  to  be  two-dimensional. 
Furthermore,  the  pressure  profiles  are  quite  similar  in  shape' showing  re¬ 
latively  little  difference  between  the  two  measuring  stations.  It  is  in¬ 
ferred  from  this  that  the  flow  is  typical  of  a  highly  developed,  turbulent 
nozzle  wall  boundary  layer. 

v  *  o 

The  static  pressure  profiles  measured  at  stations  X  “  112,  137,  and  160 
inches  from  the  throat  are  shown  in  FI  .e  11.  The  data  at  X  *  112  and 
137  inches  were  obtained  from  single  t  as  at  these  stations  while  the  X  * 

160  data  shows  the  range  of  measurements  obtained  from  six  runs.  Also  shdwn 
are  the  range  of  wall  pressure  measurements  observed  during  the  run  made,  at 
each  station.  The  wall  measurements  were  made  with  a  transducer  identical 
to  that  used  for  the  static  pressure  probe.  Although  the  difference  in 
pressure  at  successive  stations  is  only  slightly  larger  than  the  experi¬ 
mental  accuracy,  the  data  indicate  a  small  favorable  longitudinal  pressure 
gradient  of  approximately  -  0.005  mm  Hg/inch  near  the  wall.  In  the  free- 
stream  and  for  the  wall  measurements  at  the  X  *  112  and  160  inch  stations, 
the  longitudinal  variation  in  pressure  is  within  the  scatter  observed  from 
run.  As  described  in  Section  VI,  the  static  probe  measurements  are 
subject  to  large  viscous  interaction  corrections.  The  corrected  data  show 
a  pressure  gradient  normal  to  the  wall  that  is  20  times  greater  than  the 
longitudinal  gradient  indicated  above. 

A  further  indication  of  the  two- dimensionality  of  the  boundary  layer  flow, 
which  also  demonstrates  the  reproducibility  of  the  measurements,  is  shown’ 
in  Figure  12.  The  envelope  of  measurements  obtained  by  USC/JPL  showing  the 
small  spanwise  variation  of  pitot  pressure  through  the  boundary  layer  is 
presented,  together  with  the  results  of  a  pitot  pressure  survey  made  six 
months  later  during  the  present  study.  The  pitot  profile  obtained  from 
the  center  probe  is  denoted  Pt2«  The  data  indicate  that  spanwise  variations 


SECTION  VI 


FINAL  TEST  DATA 


Before  discussing  the  method  of  data  reduction  and  the  final  results  for 
the  mean  flow  profiles,  the  raw  data  is  presented  in  order  to  demonstrate 
again  the  reproducibility  of  the  measurements  and,  secondly,  to  clarify 
the  data  reduction  procedure.  The  calculation  of  the  mean  f low  properties 
is  based  on  the  measured  pitot  pressure,  static  pressure, . and  tofai^tem- 
perature.  Each  of  these  parameters  is  discussed  in  d^tai^  in  the  following 
paragraphs . 

6.1  PITOT  PRESSURE  MEASUREMENTS 

''  s  • 

Figure  13  shows  the  pitot  profiles,  denoted  as  Pt2,  obtaine^ with  the  cen¬ 
ter  pitot  probe  during'  a  total  of-yfive  runs.  Included  also  are  the  mea¬ 
surements  obtained  with  the  miniature  probe  located  at  one  of  the  outboard 
petitions.  These  points,  denoted  as  Pt3,  are  restricted  to  Y  distances 
between  0  and  1.5  inches.  As  shown  in  the  insert  cf  Figure  13,  the  dif¬ 
ferences  between  Pt2  and  Pt3  are  not  significant  anJ  two  profiles  are 
nearly  identical  fotY  >  0*5  inches.  Because  of  the  excellent  agreement 
between  the  results  jdf  different  runs,  the  data  from  all  runs  were  combined 
to  a  single  pitot  pressure  profile  for  purposes  of  data  reduction. 

Values  of  Pt3  were  used  for  0  £  Y  ^  0.8  inches  and  Pt2  data  was  used  between 
Y  "*  0.8  inches  and  the  freb-stream.  4 


6.2  TOTAL  TEMPERATURE  MEASUREMENTS 


The  measured  total  temperature  profiles  for  three  runs  are  shown  in  ‘Figure  . 
14  which,  for  comparison,  also  includes  the  data  obtained  by  the  USC/JPL 
team. (3)  Again  excellent  agreement  between  runs  is  observed  with  differ¬ 
ences  in  the  data  attributed  primarily  to  variations  in  the  supply- temper¬ 
ature.  The  insert  in,JFigure  14  shows  the  temperature  profiles'  plotted  to 
an  expanded  Y  scale  in  the  vicinity  of  the  wall  fhere  agaife,  except . for 
several  obviously  eri^>ue<Jftts  ppints  on  run  10,^ the  agreetifent  between  runs 
is  very  good.  %y‘  - 

r  • 

Similar  to  the  pitot  pressure  data,  the  data  from  alb  three  runs  was  com¬ 
bined  to  yield  aCStjaple  total  temperature  profile.  However,  since  the  Y 
position  of  the  thermocouple  did  not  coincide  with  that  of  the  miniature 
pitot  probe,  the  temperature  data  for  Y  <  0.8  inches  was  represented  by 
the  following  expressions: 


0.060  -  Y  ^  0.275  inche^ 


T  (°K)  *  303  •+  2534  Y  -  ,L$774  Y2  +  70744  Y3  -  127796  Y4  +  8896  Y5 

m  (7) 


T  (°K) 
tn 


0.275  ^  Y  -  0.800  inches 

460  +  45.7  (Y  -  0.275) 


5/ 
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POLYNOMIAL  CURVE 
FIT  TO  DATA 


FIGURE  14.  MEASURED  STAGNATION  TEMPERATURE  PROFILE  USED  FOR 
DATA  REDUCTION 


The  resulting  curve  fit  is  shown  in  Figure  14.  The  determination  of  the 
actual  total  temperature  for  Y  £  0.06,  where  temperature  measurements  were 
not  available,  is  described  ip  Section  VIII. 

6.3  STATIC  PRESSURE  MEASUREMENTS 

The  static  pressure  measurements  acquired  from  six  runs  are  plotted  as  a 
function  of  Y  in  Figure  15.  In  contrast  to  the  other  data,  the  static 
pressure  exhibited  significant  variation  on  a  given  run  as  well  as  from 
run  to  run.  As  shown  in  Figure  15,  the  maximum  deviation  in  static  pres- 
sure  for  the  six  tests  was  ±  7X,  which  was  an  order  of  magnitude  greater 
then  the  variations  in  supply  conditions.  However,  a  greater  precision 
in  the  static  pressure  measurement  should  not  be  expected  since  the  nominal 
pressure  was  less  than  IX  of  the  range  for  which  the  transducer  was  de¬ 
signed  to  operate.  Consequently,  it  was  decided  to  represent  the  measured 
static  pressure  profile  by  the  following  expression,  determined  from  a 
least  squares  fit  of  the  data: 

0  ^  Y  s  5.75  inches 

^sm^ini1  *  1.642025  -  0.197563  Y  +  0.130692  Y 2  -  0.026982  Y3  +  0.002297  Y^ 
-  0.000069  Y5 

(8) 

Y  >  5.75  inches 

p  (™n  Hg)  *  1.77 

sm 

The  resulting  curve  fit  is  shown  in  Figure  15. 

Although  the  static  pressure  measured  in  the  free-stream  was  1.77  cm  Hr 
the, static  pressure  Fgl,  calculated  from  the  free-stream  pitot  pressure 
and  the  supply  pressure,  was  approximately  1.1  cm  Hg.  The  free-stream 
Mach  numbers  determined  Irom  various  combinations  of  P  F  -  P  and  P 
are  shown  below.  '  sin*  si»  t  *o 


PRESSURE  RATIO 


psn/pt 


MACH  NUMBER 

*. 


Psra^Po  8-8 

Psi/Pt  (or  Psi/po>  9.4 

Since  the  nominal  Mach  number  corresponding  to  the  tunnel  geometry  (deter¬ 
mined  from  a  calibration"  acquired  through  long  term  use  of  the  tunnel) 
was  9.5  and  since  both  the  pitot  and  supply  pressures  are  accurate  within 

bv’p  Wa™CO™i“d*d  th3t  thC  aCtual  free-stream  static  pressure  was  given 
by  Psl.  The  difference  between  Fsm  and  Psi,  therefore,  has  been  attributed 
o  viscous  interaction  effects  on  the  static  pressure  probe.  Since  the 
interaction  effects  vanish  near  the  wall,  it  is  seen  that  a  significant 


f 


FIGURE  15.  MEASURED  STATIC  PRESSURE  PROFILE  USED  FOR  DATA  REDUCTION 


t 


pressure  difference  of  unknown  origin  existed  across  the  boundary  layer. 
Because  of  the  magnitude  of  the  pressure  difference  (approximately  50%)  it 
was  considered  necessaty  to  account  for  the  variation  in  pressure  normal 
to  the  wall  in  the  data  reduction  process. 

(7)  « 

Behrens,  using  a  static  probe  very  similar  in  geometry  to  that  used  in 
the  present  study,  developed  an  empirical  correlation  of  viscous  interaction 
effects  in  the  form 


where : 


P  /P  . 

so  si 


1  ♦  CXX  +  C2x' 


v  =  ^x2/( C/Rex)* 

C  =  (dw/^)(Tac/Tw) 


and  the  Reynolds  number  is  based  on  the  distance  between  the  tip  of  the 
probe  and  the  pressure  orifice,  and  subscripts  *  and  w  denote  free-stream 
and  adiabatic  wall  conditions  respectively.  In  our  experiments  the  maxi¬ 
mum  value  of  £  was  found  to  be  equal  to  1.05;  therefore  C  was  set  equal  to 
unity.  Behrens  data  was  restricted  to  X  <  1.6,  which  is  only  slightly 
greater  than  one-hal'  the  maximum  value  of  \  in  the  present  tests.  While 
the  results  of  Rogers,  et  al,  '  extends  to  larger  values  of  \,  their  data, 
which  was  obtained  with  a  large  cone  angle  probe,  fails  to  account  for  the 
observed  pressure  difference.  This  is  in  accordance  with  the  results  of 
Wagner  and  Watson ( who  found. that  viscous  correstions  diminish  w/th  in¬ 
creasing  cone  angle.  Consequently,  viscous  corrections  which  cover  the 
range  of  conditions  encountered  in  the  present  experiments  were  determined 
by  supplementing  Behrens  data  with  data  acquired  during  calibration  of: 

(1)  the  hot  wire  at  Mach  9.5,  and  (2)  the  total  temperature  probe  at  Mach  4. 

/ 

As  described  earlier,  calibration  was  accomplished  by  locating  the  probe 
in  the  known  free-stream  flow  at  fixed  Mach  number  and  varying  the  supply 
.pressure  to  change  the  Reynolds  number.  The  pitot  and  supply  pressures 
were  used  to  determine  the  Mach  number  and  actual  static  pressure  as  well 
as  other  flow  properties.  From  this  information  and  the  measured  static 
pressure,  P-  /Pgj  and  \  could  be  calculated.  The  results  are  shown  in  . 
Figure  16,  which  includes  also  a  plot  of  the  expression 

F  /P  .  =  1  +  0.0457  \  +  0.0622  v  (9) 

sm  si 


which  was  fit  to  the’data^ 
6.4  1700  CM  HG  TEST  DATA 


taby 


the  method  of  least  squares. 


Tbe  raw  data  obtained  at  PQ  =  *1700  cm  Hg  are  shown  in  Figure  17.-  The  data 
. .e  similar  to  that  acquired  at  3200  cm  Hg,  with,  differences  between  the 
two  attributed  to  the  increased  boundary  layer  thickness  corresponding  to 
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BEHRENS 


FIGURE  If).  VISCOUS  INTERACTION  CORRECTION  FOR  STATIC  PRESSURE  PROBE 


the  seller  Reynolds  maker  and  to  the  reduction  in  the  absolute  value  of 
the  supply  pressure.  Heasureaents  were  not  aade  for  T  <  1.0  inch  since  the 
Reynolds  nunbers  near  the  wall  are  considerably  less  than  the  range  over 
which  the  hot-wire  aneaoaeter  was  calibrated.  Use  of  the  wire  in  this 
repfon  would  have  introduced  large  uncertainties  in  interpretation  of  the 
hotwire  data. 
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SECTION  VII 


DATA  REDUCTION  PROCEDURE 


7.1  INPUT  DATA 

Calculation  of  aean  flow  properties  were  carried  out  at  each  Y  position 
for  which  a  aeasureient  of  the  pitot  pressure  was  aade.  The  input  data 
required  for  the  calculation  and  the  fora  of  the  data  are  listed  below. 


INPUT  QUAKT1TY 

SYMBOL 

FORM  OF  DATA 

Pitot  Pressure 

Pt2 

Tabular,  versus  Y 

Measured  Total  Teaperature 

T 

m 

0  <  Y  £  0.80,  Eqn.  (7) 

Y  >  0.80,  Tabular  versus  same 

Y  as  for  Pt 

Measured  Static  Pressure 

Psn 

Eqn.  (8) 

Total  Teaperature  Probe 
Recovery  Factor 

- 

Eqn.  shown  in  Figure  9 

Static  Pressure  Probe  Viscous 
Interaction  Correction 

P  /P  . 

so  si 

Eqn.  (9) 

As  mentioned  earlier,  for  Y  less  than  0.80  inches  the  Y  position  of  the 
total  teaperature  probe  did  not  coincide  with  that  of  the  pitot  probe. 
Therefore,  for  convenience  in  the  calculation  was  represented  by  Eqn.(7) 
for  Y  2  0.80.  For  Y  >  0.80  inches  the  Y  positions  of  the  two  probes  were 
identical  and  T—  was  input  to  the  calculations  in  tabular  fora. 

7.2  MEAN  FLOW  COMPUTER  PROGRAM 

Calculation  of  the  aean  flow  properties  was  programmed  in  Fortran  for  the 
Honeywell  Model  615  Coaputer  Systea.  Since'  bqth  the  measured  total  tem¬ 
perature  and  static  pressure  had  to  be  corrected  for  probe  effects,  the 
program  includes*  iterative  routines,  for  finding  the  actual  value  of  these 
parameters.  Referring  to  Figure  18,  the  pitot  pressure  and  P^  are  used 
to  find  the  Mach  nunber  M  and  the  latter  two  then  combined  to  calculate 
the  local  stagnation  pressure,  PQ.  The  Mach  number  is  also  used  with  the 
measured  total  teaperature  to  find  the  static  streae  temperature,  T  .  With 
the  teaperatures  known  the  corresponding  viscosities  are  evaluated,  and  the. 
Reynolds  numbers  based  on  the  appropriate  temperature  and  the  characteristic 
prob*  diaension  can  be  calculated.  The  Reynolds  nuober  for  the  static 
pressure  probe  is  combined  with  the  Mach  number  to  find  the  viscous  inter¬ 
action  paraaeter  \  and,  by  means  of  Eqn.  (9),  the  first  order  approximation 
to  the  static  pressure  Ps..  The  procedure  is  then  repeated  n  times  until 


REPEAT  ITERATION 


FIGURE  18.  FLOW  CHART  FOR  MEAN  FLOW  COMPUTER  PROGRAM 
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w  i 


the  absolute  value  of  t  1)  -  PS£(n)  satisfies  the  convergence  criteria 

€.  Similarly  the  total  temperature  probe  recovery  factor,  T),  is  found  from 
the  probe  Reynolds  number,  yielding  the  first  order  approximation  to  the 
local  total  temperature, » T0.  This  calculation  is  also  repeated  until* , 

|T0(n  +  1)  -  Tc(n)|  <  e.  in  both  eases  the  convergence  parameter  e  was  0.1Z. 

Results  of  the.iterW^ve  calculation  to  determine  P  at  Y  *  6.016  inches  are 
presented  in  Figure  19 T  which  indicates  that;  12  iterations  are  required  to 
satisfy  the  desired  convergence.  Since  the  viscous  correction  is  greatest 
owr  the  free-stream,  this  example  is  typical  of  the  maxiaum  number  of 
iterations  required.  Convergence  of  the*  total  temperature  generally  occur¬ 
red  more  rapidly,  although  convergence  pf  both  pari-*eters  was  required  be-  >“ 
fore  cobtiouing  the  calculation  procedure.-  , 

‘  .  ■ 

Once  T0  and  -PS£  were  determined  the  remaining  boundary  layer  properties 
were  calculated,  i;e.,  the  velocity  u,  the  density  P,  etc.,  using  standard 
'formulae  for  compressible  flow.  In  addition  to  the  local  properties, 
several  characteristic  boundary  layer  thicknesses,  including  the  displace¬ 
ment  thickness  5*,  and  the  momentum  thickness,* s,  were  determined  frors 
numerical  integration  across  the  boundary  layer.., 


t 


_  % 

/*• 


a;  ' 


r 


■  > 


i 


/- 


U 


-36- 


Q 


7 


SECTIOX  VIII 
RESULTS 

•+-~*su  *  j. 

«  ’ 

Hie  resists  of  the  mean  flow  calculations  at  PQ  =  3200  cm  Hp  are  suimaarized 
in  Table  I  which  list?  the  test  conditions,  the  characteristic  boundary 
layers  thicknesses  described  in  Paragraph  7.2,  and  a  tabulation  of  the 
oean  flow  profiles.  The  conventional  boundary  layer  thickness  ',  which  is 
defined  as  the  distance  from  the  wall  where  the  velocity  '.s  99%  of  the 
free-stream  velocity,  has  been  included  although  in  the  present  case  its 
significance  is  somewhat  vague.  This  arises  from  the  fact  that  relatively 
few  data  points  were  obtained  for  Y  >  5.0  inches,  thereby  precluding  a  pre¬ 
cise  determination  of  the  shape  of  the  velocity  versus  Y  curve  near  the 
edge  of  the  boundary  layer.  For  this  reason  Y,  rather  than  Y/-,  has  been 
used  in  subsequent  presentation  of  the  results,  even  though  the  results 
have  been  nondimensionalized  using  the  values  of  the  flow  properties  at 
Y  *  i  as  reference  paraaeters.  It  should  be  pointed  out  that  this  repre¬ 
sents  an  arbitrary,,  although  consistent,  selection  of  reference  parameters 
since  the  velocity  and  thermal  boundary  layers  are  of  different  extept, 
with  the  latter  being  slightly  larger.  Furthermore,  because  of  the  pres¬ 
sure  gradient  normal  to  the  wall,  the  fluid  density  actually  overshoots 
its  free-streaa  value  within  the  boundary  layer. 

8.1  *AK  FLOW  PROFILES 

'  ■  \  — 

8.1.1  STATIC  PRESSURE  PROFILE 

(9) 

The  calculated  static  pressure  profile  is  shown  in  Figure  13.  Fiore 
and  Fischer  et  al'*®'  have  also  observed  pressure  gradients  normal  to  the 
wall  in  hypersonic  boundary  layers,  the  former  at  Mach  12  in  air  and  with 
the  Reyrolds ^number  based  on  momentum  thickness,.  Re^,  on  the  order  of  1000 
and  the  latter  at  Mach  21  in  helium  with  Re^  ~  10,000.  In  Fiore's  casej^) 
the  pressure  gradient  was  smaller  than  observed  here,  with  the  pressure 
difference  equal  to  20%  of  the  free-stream  pressure,  but  was  still  larger 
thajn  the  longitudinal  pressure  gradient  neasurecl  in  his  tunnel  and  was  (10) 

found  to  increase  with  increasing  Re^.  In  the  experiments  of  Fischer  et  al 
the  static  pressure  at  the  wall  was  40%  greater  than  the  free  stream  pres¬ 
sure,  similar  to  our  results.  Fisher  et  al^10*  also  examined  the  results 
of  a  number  of  other  investigators  who  reported  an  increase  in  measured 
wall  pressure  over  the  corresponding  free-stream  value  for  several  gases 
and  variety  of  aerodynamic  bodies,  including  cones,  wedges,  and  tunnel 
Vails.  They  found  a  consistent  trend  of  increasing  wall- to- free- stream 
static  pressure  ratio  with  increasing  Xach  number  which  was  tentatively 
attributedto  turbulent  flow  effects.  \ 

t  ,  s 

8.1.2  TEMPERATURE  PROFILES 

Figure  20  shows  the  ratios  of  local  stagnation  and  static  temperature  to 
the  stagnation  temperature  at  the  edge  of  the  boundary  layer,  1^,  (i.e., 
its  value  at  Y  *  j)  plotted  against  Y.  The  insert  in  Figure  20  shows  the 
region  near  the  wall  plotted  to  an  expanded  Y  scale.  Since  temperature 
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TABLE  I 


SUMMARY  OF  RESULTS 


p 

o 

3200  cm  Hg 

P 

se 

1.15  on  Hg 

Re- 

36,800 

7 

792°K 

u 

1225  cs/sec 

k 

6.01 

inches 

oe 

e 

-§r 

T 

V 

304  1C 

M 

e 

9.37 

6 

2.29 

•o 

inches 

T 

e 

43°K 

Re 

o 

127,080 

c 

.29  : 

Inches 

No. 

v (inches) 

M 

u/u 

T  /T 

T/T 

P/P 

e 

o  oe 

oe 

e 

11 

0.010 

0.541 

0.161 

0.448 

0.424 

0.180 

2 

0.015 

0.541 

0.166 

0.477 

0.452 

0.170 

3 

0.015 

6.663 

0.201 

0.477 

0.438 

0.175 

4 

0.018 

0.760 

0.236 

0.493 

0.457 

0.168 

5 

0.022 

0.607 

0.189 

0.510 

0.475 

0.161 

6 

0.031 

0.842 

0.268 

0.545 

0.477 

0.160 

7 

0.037 

0.843 

0.276 

0.564 

0.504 

0.151 

8 

0.049 

1.130 

0.356 

0.591 

0.471 

0.163 

9 

0.054 

1.083 

0.347 

0.600 

0.486 

0.157 

10 

0.057 

1.354 

0.410 

0.605 

0.442 

0.176 

11 

0.063 

1.131 

0.363 

0.613 

0.488 

0.156 

•  12 

0.070 

1.260 

0.396 

0.616 

0.468 

0.163 

13 

0.075 

1.261 

0.398 

0.622 

0.472 

0.161 

14 

0.088 

1.497 

0.452 

0.625 

0.431 

0.176 

15 

0.110 

1.670 

0.490 

0.633 

0.406 

*0.186 

16 

0.122 

1.672 

0.492 

0.639 

0.410 

0.184 

17 

0.172 

1.912 

0.537 

0.645 

0.373 

0.201 

18 

0.209 

2.004 

0.553 

0.649 

0.360 

0.207 

19 

0.219 

2.005 

0.554 

0.652 

0.361 

0.206 

20 

0.222 

2.055 

0.561 

0.650 

0.353 

0.211 

21 

0.273 

2.155 

0.579 

0.660 

0.342 

0.216 

22 

0.307 

2.193 

0.585 

0.661 

0.337 

0.219 

23 

0.322 

2.229 

0.590 

0.661 

0.331 

0.222 

24 

0.324 

2.206 

0.587 

0.662 

0.335 

0.220 

25 

0.370 

2.311 

0.601 

0.661 

0.320 

0.229 

26 

0.407 

2.347 

0.606 

0.663 

0.315 

0.232 

27 

0.422 

2.370 

0.609 

0.663 

0.312 

0.234 

28 

0.472 

2.438 

0.618 

0.665 

0.304 

0.239 

29 

0.521 

2.463 

0.622 

0.668 

0.302 

0.240 

30 

0.624 

2.592 

0.6.18 

0.671 

0.286 

0.251 

31 

0.722 

2.696 

0.650 

0.675 

0.275 

0.260 

32 

0.808 

2.869 

0.666 

0.675 

0.255 

0.279 

33 

0.815 

2.888 

0.671 

0.681 

0.255 

0.279 

34 

0.913 

2.999 

0.688 

0.696 

0.249 

0.285 

35 

1.064 

3.206 

0.701 

0.690 

0.225 

0.312 

36 

1.362 

3.584 

0.734 

0.709 

0.199 

0.353 

37 

1.515 

3.743 

0.753 

0.727 

0.191 

0.366 

0 
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No. 

V (inches) 

M 

38 

«■  (.662 

3.975 

39 

1.962 

4.348 

40 

2.056 

4.459 

41 

2.114 

4.587 

42 

2-268 

4.882 

43 

2.562 

5.353 

44 

2.712 

5.672 

45 

2.761 

5.728 

46 

2.864 

6.000 

47 

3.012 

6.267 

48 

3-116 

6.731 

49  / 

3.308 

-  6.939 

50 

3.311 

6.901 

3.508 

7.431 

52 

3.612 

7.755 

53 

3.708 

7.854 

54  v 

3.712 

7.974 

55 

3.811 

8.226 

56 

3.910 

8.277 

57 

3.911 

8.477 

58 

4.013 

8.576 

59 

4.015 

8.591 

60 

4.107 

8.567 

61 

4.114 

8.746 

62 

4.215 

8.833 

63 

4.309 

8.800 

64 

4.319 

8.903 

65 

4.416 

8.993 

66 

4.509 

8.922 

67 

4.510 

9.008 

68 

4.812 

9.157 

69  ' 

5.H1 

9.202 

70 

5.402 

9.294 

71 

6.016 

9.373 

72 

7.015 

9.345 

73 

8.016 

9.364 

t 


TABLE  1  -  (Continued) 

u/u^ 

*T  /T 
o  oe 

T/T 

oe 

0.765 

0.729 

0.175 

0.792 

0.751 

0.157 

0.797 

0.752 

0.151 

0.811 

0.769 

0.148 

0.823 

G.776 

0.134 

0.848 

0.799 

0.119 

0.863 

0.815 

0.110 

0.860 

0.807 

0.107 

0.874 

0.823 

0.100 

0.884 

0.834 

0.094 

0.898 

0.847 

0.084 

0.904 

0.853 

0.080 

0.909 

0.863 

0.082 

0.920 

0.874 

0.073 

0.927 

0.881 

0.068 

0.929 

0.882 

0.066 

0.934 

0.891 

0.065 

0.939 

0.896 

0.062 

0.942 

0.902 

0.061 

0.943 

0.901 

0.059 

0.948 

0.908 

0.058 

0.951 

v  6.915 

0.058 

0.946 

0.904 

0.058 

0.952 

0.914 

0.056 

..  0.958 

0.923 

0.056 

0.957 

0.923 

0.056 

0.960 

0.927 

0.055 

0.964 

0.933 

0.054 

0.964 

0.935 

0.055 

0.964 

0.933 

0.054 

0.976 

0.954 

0.054 

0.984 

0.970 

0.054 

0.990 

0.981 

0.054 

1.001 

1.001 

0.054 

1.012 

1.025 

0.055 

1.010 

1.020 

0.055 

c/P 


0.398 

0.444 

0.462 

0.471 

0.316 

0.581 

0.625 

0.643 

C.677 

0.716 

0.785 

0.820 

0.802 

0.883 

0.929 

0.948 

0.954 

0.988 

0.992 

1.023 

1.0*9 

1.020 

1.038 

1.046 

1.046 

1.045 

1.053 

1.056 

1.048 

1.061 

1.050 

1.032 

1.027 

0.999 

0.958 

0.966 
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measurements  could  not  be  made  for  Y  <  0 . 062)  inches,  the  variation  of  total 
temperature  in  this  region  was  determined  by  fitting  a  third  degree  poly- 
nominal  curve  to  the  data  points  at  Y  *  0  (where  Tq/T^  ”  Tw/Toe  )  and  Y 
*  0.063  through  0.122.  Sio^e  the  viscous  correction  to  Psm  was  zero  near 
the  wall,  Psi  could  be  determined  from  Eqn.  (8)  and  the  remaining  flow 
properties  evaluated  explicitly.  The  resulting  Ts/Toe  profile  shows  the 
anticipated  static  temperature  overshoot  and  .udicates  that  the  maximum 
static  temperature  occurs  at  Y  as  0.04  inches  wiser  it  exceeds  the  wall 
temperature  by  approximately  80°K. 

8.1.3  MACH  NUMBER  PROFILE 

rf 

The  Mach  number  profile  is  shown  in  Figure  21  which  indicates  a  gradual 
decrease  in  Mach  numbe**  from  its  free-stream  vsiue  to  M  *  9  at  Y  =  4.3 
inches.  Between  Y  ■  4.5  and  Y  “  2.5  inches,  the  Mach  number  decreases 
rapidly  in  a  nearly  linear  fashion  to  M  *  5,  followed  by  a  slower,  but 
again  nearly  constant  rate  of  decrease  to  M  *  2  at  Y  =  0.20  inches.  The 
sonic  line  occurs  at  approximately  0.05  inches,  so  that  subsonic  flow  is 
restricted  to  a  very  narrow  layer  adjacent  to  the  wall. 

8.1.4  DENSITY  PROFILE 

The  density  profile,  shown  in  Figure  22,  indicates  that  at  Y  *  4.5  inches 
the  density  overshoots  its  free-stream  value  by  101.  This  reflects  the 
fact  that  as  the  wall  is  approached,  the  static  pressure  begins  to  increase 
over  its  free-stream  value  at  Y  *  6.0  inches,  while  the  static  temperature 
does  not  increase  appreciably  until  Y  -  4.0  inches.  The  minimum  density 
corresponds  to  the  static  temperature  overshoot  and  is  about  25%  less  than 
its  value  at  the  wall.  ' 

8.1.5  VELOCITY  PROFILE  1 

The  velocity  profile  is  presented  in  Figure  23  which  shows  that  over  one- 
half  the  velocity  change  across  the  boundary  layer  occurs  in  the  imediate 
neighborhood  of  the  wall,  i-.e.,  u/ue  decreases  from  0.55  to  zero  between 

Y  *  0.2  inches  and  Y  *  0.  This  region  has  been  plotted  with  an  enlarged 

Y  scale  in  the  insert  in  Figure  23,  where  it  is  noted  that  there  is  no 
evidence  of  the  linear  variation  of  velocity  generally  attributed  to  the 
viscous  sub-layer.  This  subject  is  discussed  further  in  Section  IX. 

8-2  COMPARISON  WITH  HOT-WIRE  DATA 

The  availability  of  the  hot-wire  data  provides  a  means  for  partially  check¬ 
ing  the  mean  flow  measurements.  At  any  point  in  the  boundary  layer,  the 
recovery  temperature,  Tw,  of  the  wire  and  the  Nusselt  number  can  be  de¬ 
termined  from  t’je  overheat  traverse  made  during  the  turbulence  measure¬ 
ments  and,  using  the  Hu  versus  ReQ  calibration  data  shown  in  Figure  8, 

Re0  can  be  found.  Then,  with  ReQ  and  known,  Tjcan  be  found  from  Figure 
9  and  the  cotal  temperature  evaluated.  The  calculation  procedure  is  des¬ 
cribed  in  Appendix  A. 

i1 


Figure  24  shorn  a  comparison  of  T0  profiles  obtained  from  the  hot-wire  and 
thermocouple  measurements.  For  wire  A13-1,  the  agreement  between  the  two 
measurements  is  excellent,  with  the  hot-wire  results  generally  no  more  than, 
5-10°C  below  the  thermocouple  data .  A  larger  discrepancy  is  observed  for 
wire  6-3,  although  the  agreement  with  the  thermocouple  results  is  still 
encouraging .  In  this  case,  the  hot  wire'  results  are  about  20°C  higher  than 
those  obtained  from  the  thermocouple  measurements  in  the  outer  half  of  the  i 
boundary  layer,  with  the  discrepancy  increasfftg  to  50°C  between  Y  *  1  and 
2  inches,  then  gradually  4isioi shing  as  the  wall  is  approached. 

Figure  25  presents  a  comparison  of  the  profiles  of  the  unit  Reynolds  num¬ 
ber  per  inch  based  on  stagnation  temperature,  Re0'-  Although  the  shape  cf 
the  curves  are  similar,  ^ -Q'  obtained  with  wire  A13-1  exceed  those  deduced 
from  the  mean  flow  measurements  near  the  wall,  while  in  the  outer  half  of 
the  boundary  layer  the  converse  is  true.  In  spite  of  the  large  differences 
between  the  two  results,  which  reach  a  maxi  nun  of  501,  the  agreenent  is 
again  encouraging.  This  arises  from  the  fact  that  snail  errors  in  XuQ  lead 
to  errors  in  ***>■  which  are  2  to  3  times  larger.  This  is  indicated  in 
Figure  8  which  shows  also  the  typical  scatter  to  be  expected  in  the  deter¬ 
mination  of  Xu0. 

*  <• 

The  Re0*  profile  obtained  with  wire  6-3  shows  taich  larger  deviations  from 
that  deduced  from  the  mean  flow  measurements,  with  differences  as  large  as 
a  factor  of  ten  observed  near  the  wall.  This  prompted  a  re-examination  of 
the  Mu  and  "  versus  Re  calibrations  froa  which  it  was  concluded  that  the 
data  points  at  ReQ  *  0.026  and  0,11  (see  Figure  8)  may  be  in  error.  Con¬ 
sequently  ,  new  correlations  were  curve  fit  to  the  remaining  five  data  points 
and  the  hot-wire  data  reduction  was  repeated.  The  resulting  Re0*  profile. 
Figure  25  shows  auch  better  agreement  with  the  results  obtained  from  both 
the  mean  flow  measurements  and  wire  A13-1-  In  addition,  the  Tq  profile 
obtained  with  the  modified  hot-wire  calibration  curves,  presented  in  Figure 
17,  indicates  slightly  better  agreement  with  the  thermocouple  measurements, 
particularly  near  the  wall  and  in  the  region  from  Y  *  1  to  2  inches  where 
the  difference  between  the  two  'sets  of  results  is  reduced  from  50  to  30-40°C. 
• 

8.3  TEST  RESULTS  AT  PQ  -  1700  CM  BG 

The  calculated  static  pressure  profile  is  shown  ip  Figure  17  where  it  can 
be  compared  to  the  measurements.  The  pressure  variation  is  similar  to  that 
observed  at  FQ  *  3200  cm  Eg  although,  because  of  the  reduced  Reynolds  num¬ 
ber,  the  viscous  corrections  at  the  outer  edge  of  the  boundary  layer  are 
slightly  larger. 

The  static  and  stagnation  temperature  profiles  and  the  velocity  and  Mach, 
number  profiles  are  shown  in  Figures  26  and  27,  respectively,  where,  for 
comparison,  the- results  at  Fa  ”  3200  cn  Eg  have  been  included,  ftie  to  the 
lack  of  data  near'  the  edge  of  the  boundary  layer,  an  accurate  calculation 
of  5  could  not  be  made.  However,  because  of  the  similarity  in  the  boundary 
layer,  profiles  and  the  known  equivalence  of  the  free-str ears  conditions,  it 
was  assumed  that  the  values  of  ue  and  T^  at  the  edge  of  th^  boundary  layer 
were  identical  to  those  at  P„  •  3200  cn  Eg  and  these  values  were  used  to 
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non- dimens localize  the  PQ  -  1700  ca  Hg  data.  The  resulting  non-dimensional 
profiles  are  in  good  agreement  with  those  obtained  at  the  higher  pressures 
and  the  differences  between  the  two  could  be  reduced  by  normalizing  Y  with 
respect  to  a  characteristic  boundary  layer  thickness. 

The  lack  of  data  also  prevented  a  precise  calculation  of  Re& .  To  a  first 
approximation,  however,  this  parameter  should  scale  as  the  ratio  of  the 
pressures  so  that  for  the  PQ  »1700  cm  Hg  tests.  Re*  as  20,000; 
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SECTION  IX 


a 

DISCUSSION  OF  RESULTS 


9.1  COMPARISON  WITH  OTHER  EXPERIMENTS 


Figure  28  shows  the  local  static  ^esqjerature  ratio  within  the  boundary 
layer,  T/Te,  plotted  against  the  loS'ai' velocity  ratio,  u/u£.  For  compari¬ 
son  the  Crocco  relation  « 


T/T 

e 


T  /T  +  (1  +  (v-1) 
w  e 


M  2/2  -  T  /T  )  u/u 
ewe  e 


0-1) 


(u/u  )2/2 

f" 


(10) 


which  is  valid  for  arbitrary  heat  transfer  and  zero  lohgitudinal  pressure 
gradient,  is  included.  While  there  is  good  agreement  between  experiment 
and  theory  in  the  sub-layef  region,  the  data  lies  significantly  below  the 
Crocco  relation  in  the  outer  portion  of  the  boundary  layer.  This  _is  illus¬ 
trated  more  clearly  in  Figure  29,  where  the  total  enthalpy  ratio,  T  =  TQ  - 
Tw/Tpe  -  Tu,  is  plotted  versus  u/ue.  In  these  coordinates  the  Crocco  rela¬ 
tionship  is  represented  by  the  linear  function  T  =  u/ue.  Again,  in  the  sub¬ 
layer  the  data  agrees  with  the  theory,  while  in  the  outer  portion  %f  the 

boundary  layer,  corresponding  to  the_law-of-th^-wake  region,  the  data  iol- _ 

lows  closely  the  quadratic  relation  T  =  (u/u.)  .  The  interior  of  the  boun¬ 
dary  layer  -  the  so-called  iaw-of-the-wail  region  -  is  characterized  by  a 
well  defined  transition  between  the  linear  and  quadratic  behavior  of  the 
T  versus  u/u.  relationship.  *  .. 

6  1  r 

The  total  enthalpy  profile  of  Figure  29  is  typical  of  that  obtained  from 
measurements  on  nozzle  walls.  Bertram  and  Seal'  '  and  later  Bushel 1  et  al 
examined  the  results  of  numerous  experiments  obtained  for  «a  variety  of  oper¬ 
ating  conditions  including  Mach  numbers  varying  from  3  to  20,  Re^-  ranging 
from  1000  to  50,000  and  Tw/Toe  varying  from  0.2  to  1.0.  They  found  that 
ddta  acquired  from  flat  plate  measurements  generally  followed  the  linear 
Crocco  relation  while  the  nozzle  wail  data  approximated  the  quadratic  lav. 

The  behavior  of  the  nozzle  wall  data  has  been  attributed  by  those  and  other 
investigators  to  the  effects  of  the  rapid  expansion  in  the  wincT  tunnel  nozzle. 
They  suggest  that,  as  a  consequence  of  the  flow  acceleration  near  the  nozzle 
throat,  the  boundary  layer  profile  distorts  from  the  Crocco  relationship  and 
remains  out  of  equilibrium  with  the  local  edge  conditions  for  considerable 
distances  downstream.  For  hypersonic  nozzles,  however,  the  gradients  in 
the  flow  properties  are  concentrated  immediately  downstream  of  the  throat 
and  it  seems  unlikely  that  effects  of  the  past  history  of  the  flow*  would 
persist  for  very  large  distances. t*  Figure  30  shows  a  plot  of  T  versus  u/uf 
obtained  by  the  USC/JPL  investigators*3)  at  X  =  160  inches  from  the  nozzle 
throat.^  The  rg^ults  at  the  two  stations,  which  are  separated  by  approximately 


it r 


The  shift  in  the  transition  region  from  u/u  e  =  0.7  to  u/ue  =  0.3  is  due 
primarily  to  an  alternate  definition  of  the  total  temperature  probe  recovery- 
factor  used  by  USC/JPL  in  their  data  reduction.  The  data  reduction  process 
at  both  stations  was  consistent,  however,  and  the  conclusions  inferred  from 
Figure  30  £onc&rning  the  streamwise  variation  in  flow  profiles  and  are  in¬ 
dependent  of  the  probe  calibration.  . 
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FIGURE  29.  NON-DIMENSIONAL  TOTAL  TEMPERATURE  -  VTLOCITY  PROFILE 
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eight  boundary  layer  thicknesses,  are  identical  within  the  experimental  ac¬ 
curacy  and  show  no  evidence  of  relaxation  from  the  quadratic  to  the  linear 
temperature  distribution.  The  relaxation  process  was  investigated  by  v- 

Bushnell,  et  al,  who  performed  a  numerical  experiment  to  examine  the 
effects  (if  the  nozzle  acceleration  on  the  boundary  layer  structure.  Using 
nozzle  contours  for  Mach  6,  8  and  19  wind  tunnels  for  which  wall  measure¬ 
ments  are  available,  they  assumed  a  Crocco  temperature  distribution  at  the 
£hroat  and  solved  the  boundary  layer  equations  to  determine  the  resulting 
T  -  u/ue  profiles  at  downstream  stations.  Their  results  indicated  a  shift 
from  the  linear  relation  toward  the  quadratic  law,  particularly  near  the 
boundary  layer  eA'gt^,  with  increasing  X  and  then  a  return  to  the  linear  pro- 
file^  but  showed  ‘no  evidence  of  the  transition  region  observed  in  Figure  29. 
Furthermore,  their  calculations  predicted  a  more  rapid  relaxation  Lo  the 
linear  distribution  than  observed  experimentally,  although  at  Mach  6,  they 
found  that  a  distance  equivalent  to  60  boundary  layer  thicknesses  may  be 
required  for  relaxation  to  occur.  Thijs  may  explain  the  lack  of  streamwise 
variation  in  the  total  temperature  profile  indicated  in  Figure  30. 

Several  additional  comments  concerning  the  results  of  Figure  29  are  in 
order.  First,  Fiore^s1  '  measurements,  whith  were  obtained  at  Mach  12  and 
Re-.  1000,  show  a  quadratic  temperature  distribution  in  the  outer  portion 

of  the  boundary  layer.  However,  in  the  sub- layer  his  result*  can  he  ev- _ 

pressed  in  the  form  T  =  j (u/ue) ,  where  y  ~  0.7,  and  do  not  indicate  the 
transition  region  observed  in  Figure  29.  In  addition,  Fiore’s  data  shows 
no  evidence  of  the  law-of- the-w^ke  component  of  the  boundary  layer  which, 
because  of  the  small  Re-,  indicates  that  his  boundary  layer  may  not,  have 
been  fully  developed.  Second,  the  results  of  Lee,  et  al/13)  who  niade 
nozzle  wall  measur -tsonts  at  Mach  3  and  Re„  ranging  from  3000  to  30,000  are 
in  accord  with  the  general  features  of  Figure  29.  They  also  observed  a 
transition  from  the  quadratic  temperature  distribution  in  the  outer  portion 
of  the  boundary  layer  to  the  Crocco  relation  in  the  sub- layer.  Their  re¬ 
sults  indicated  further  that  at  a  given  station  the  transition  region  shifts 
closer  to  the  wall  with  increasing  heat  transfer  and  Re_ .  While  measure¬ 
ments  were  performed  at  stations  ranging  from  48  to  94  inches  from  the 
nozzle  throat  (the  inter-station  distance  is  approximately  16  boundary  layer 
thicknesses),  again  no  evidence  of  streairwise  relaxation  was  apparent. 

>1 

It  is  clear  from  the  preceding  discussion  that  the  present  result  is  typi¬ 
cal  of  those  obtained  measurements  in  wind  tunnel  wall  boundary  layers. 

The  observed  deviation  from  the  Crocco  total  temperature  distribution  im¬ 
plies  that  the  boundary  layer  is  not  yet  in  equilibrium  with  the  local  edge 
conditions.  However,  '.he  reason  for  the  departure  from  equilibrium  has  not 
yet  been  resolved. 

9.2  CORRELATION  OF  VELOCITY  PROFILE 

To  further  demonstrate  the  validity  of  the  measurements,  the  resulting 
velocity  profile,  Figure  23,  was  correlated  with  the  "universal"  incom¬ 
pressible  profile  which  is  conventionally  represented  in  the  form: 

sub- layer: 

+  + 
u  «  Y 
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and  A  and  B  are  constants  characterizing  the  velocity  profile,  C  is  a  con¬ 
stant  representing  the  strength  of  the  wake  cqmponent  of  the  boundary  layer, 
and  W  is  Cole’s^  '  tabulated  wake  function.  All  of  the  variables  in  the' 
above  relations  refer  to  tjie  incompressible  flow  f-ield  with  which  the  com¬ 
pressible  boundary  layer  is  being  compared.  The  parameters  u*  and  Y*  de¬ 
note  the  compressible  flow  properties  after  transformation  or  "stretching" 
to  account  for  the  density  variations  within  the  boundary  layer  (actually 

only  one  of  thgse  pa rameters  is  transformed^  the  other  retains-fCs  physical - 

identity).  The  friction  velocity,  u-,  is  found  either  directly  from  the 
slope  of  the  physical  velocity  profile  at  the  wall,  or  as  a  consequence  of 
the  correlation  procedure  when  it  serves  as  an  "adjustable  constant"  whose 
value  is  -manipulated  to  provide  the  best  fit  of  the  data  to  the  theory. 

In  the  case  of  a  cold  wall,  Walz(15)  and  Lee  et'T[il3)  caution  against 
using  the  measured  velocity  profile  to  determine  u-  because  of  -the  temper¬ 
ature  overshoot  which  occurs  adjacent  tothe  surface.  They  reasoned  that 
if  the  shear  stress  is  to  remain  constant  in  the  sub-layer,  then  the  velo¬ 
city  gradient  must  compensate  for  variations  in  .the  viscosity  arising  from 
the  temperature  distribution.  This  would  lead  to  a  nonlinear  velocity  pro¬ 
file  near  the  wall.  Since  velocity  measurements  very  close  to  the  surface 
are  frequently  lacking  or  unobtainable,  the  actual  velocity  profile  may  be 
obscured  by  the  missing  information.  This,  in  turn,  could  introduce  large 
errors  in  determining  du/dY. 

In  this  Context.  Walz  examined  the  cold  wall  data  of  Lobb  et  al  and 

Winkler  and  Cha  1)  who  had  fit  a  straight  line  to  their  velocity  measure¬ 
ments  in  order  to  determine  the  velocity  gradient  at  the  wall.  Comparing 
their  results  to  values  of  du/dY  calculated  from  a  similarity  law  he  had 
developed  for  evaluating  skin  friction  coefficients,  Walz  concluded  that 
the  experimentally  determined  slopes  were  too  low,  .in  some  cases  by  as 
much  as  507..  The  effect  of  cooling  on  the  sub-layer  velocity  profile  was 
investigated  further  by  Lee  et  al  (13)  who  used  the  method  of  TetervindS) 
to  calculate  the  boundary  layer  structure  for  Mach  10  and  several  values 
of  Tw/Taw,  where  subscript  aw  denotes  the  adiabatic  wall  condition.  Their 
results  showed  that  the  curved  portion  of  the  velocity  profile  extended 
closer  to  the  wall  as  Tw/Taw  decreased,  although  very  low  wall  temperatures 
are  required  for  the  effect  to  become  apparent. 


law-of- the-wa 1 1 ; 
law-of-the-wake ; 


u+  -  A  ln\  Y+ 


V 

n  B  Y 


u  -(ue/u_)  *  A  In,  (y/6)  -,  C  [2 


-  W(y/(). 


where 


+  */ 
u  -  u  /u_ 

+  * 

Y  -  Y  u  Jfv 


u_  ”  C  /0  ) 
w  w 


-57- 


l;  XZiTiZnu  acquirhed  d“rins  ,he  »— 
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temperature  overshoot  of  th*  ua n  *  hardly  conclusive,  the  static 

tor  ted  velocity  Profile  Jther  T  ‘V  P®«ible  cause  for  the  dis- 

Pitot  pressure  trlnfa  factors,  such  as  poor  resolution  of  the 

were  cm^Uy  as^elT  "f*  a"d  “»«rt.i„tieS  in  the  V  position, 

in  spite  of  the  lov  resolution  nf^h  ^  inconclusiv<!  "suit's.  For  example, 

ibie  senrally  conai“a"' 

*e  wicnin  1U7.  on  successive  runs.  FinaHv  it  ,  a  u  , 

2  :izli<  zzi, 

to  determine  du/dV  from  the  .1™.  curvTfU 

«*«  ZJzzsjts  zzi%\  - 

than  0.10  i nches  TtiTck  orabout  1  -  27  of  the  total  k  s^blayer  «  loss 

ness.  This  ls  attributed  prlmtt  ly  to  the  Urge  Re  ITT*  *" 

In  this  respect  it  in  of  i„t  •  1  rge  Re:;  °{  the  experiment. 

layer  thickness  of  «  2^"  ^  iTliT  l^  V  3  Sub' 

Ree  of?his  experiment  was  10*  -  in  r-  2  ltb°ugh«  as  indicated  earlier,  the 

study.  Kemp  and  0wens<19>  also  repo^ed^Lrsu^la^31^  °f  ^  PreSent 

«  ^  te-rirsr? 

-  -  =£ - 

the  boundary  layer.  region  in  the  outer  portion  of 

Correlation  of  the  velocity  profile  was  carrier!  r.. 

veloped  by  Maise  and  McDonalH^o)  ?  ,  d  out  “sing  a  procedure  de-  , 

by  Cran(2l>  to  include^ ^  a  f?'  adlabatic  flat  P^tes  and  extended 
fitting  the  iat.  to  th.  Cold.”al  """“on.  The  procedure  involves 
ing  Coles ^22)  has  been  “h‘Ch-  f°U°“- 


* 

u  /u. 


/  • 

2.43  In  (Y  u_/vw)  (  5.0  +  2.43  ~  W(Y/6) 
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for^f^  If  ^  generalized  velocity  proposed  by  Van  Driest*23* 
for  the  effects  of  compressibility  and  heat  transfer,  i.e.. 
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It  should  be  pointed  out  that  the  Van  Driest  formulation  assumes  that  the 
static  pressure  is  constant  in  the  boundary  layer  and  that  the  Crocco  rela- 
' icn,  Eqn.  (10),  is  valid,  tt.ro  conditions  which  are  not'met  in  the  present 
experiment.  This  must  be  recognized  when  interpreting  the  results  of  the 
correlation. 

,  Equation  (12)  includes  three  unknown  constants,  u_ ,  5,  and  ",  a  parameter' 
N"'^ep  resen  ting  the  strength  of  the  wake  component  of  the  boundary  layer.  By 
applying  Eqn.  (12)  at  the  edge  of  the  boundary  layer,  "  can  be  expressed 
in  terms  of  £  and  u-  and  the  edge  conditions.  This  reduces  the  number  of 
unknown  constants  in  Eqn.  (12)  to  two,  whose  values  are  adjusted  until  the 
data  fits  the  equation  such  that  the  rras  error  is  a  minimum.  Data  points 
for  Y+  *  50  and  Y  >  0.9  :  are  excluded  since  in  the  former  caise  the  experi¬ 
mental  errors  are  usually  large  while  in  the  latter  case  the  theory  deviates 
from  the  data.  For  convenience  in  the  calculations  Cole's  wake  function 
W  is  expressed  as 

M  -  2  sin2  (-Y/25)  (13) 


Af*er  at  fit  of  the  data  is  achieved,  **  and  cf,  the  skin  friction  coefficient, 
are  determined  from.i  and  u  . 

The  results  of  the  correlations  are  shewn  in  Figure  31.  The  value  of  “  was 
1.4,  slightly  more  than  twice  the  value  of  0.60  normally  found  in  correlat¬ 
ing  zero  pressure  gradient  flat  plate  data  and  typical  of  flows  with  strong 
adverse  pressure  gradients . On  thf  other  hand,  the  skin  friction  coef¬ 
ficient  was  4.6  x  10-4, .which  is  within  several  percent  of  the  theoretical 
value  determined  from  the  Karman-Schoenberr  incompressible  formula  using 
th«  Re^  listed  in  Table  I  together  with  the  Van  Driest  transformation  func¬ 
tions  and  following  th^  procedure  outlined  by  Hopkins  and  Inouye.^^  Since 
the  assumptions  of  the  Van  Driest  transformation  imply  a  specific  density 
profile,  a  fictitious  Re^  can  be  calculated  from  the  result*  of  the  velocity 
^  correlation.  The  Re..  found  in  this  manner  was  19,100,  corresponding  to  a 
101  increase  in  the  theoretical  Cf,  which  is  approximately  one-half  of  the 
•  actual  Re^  determined  from  integration  of  the  mean  flow  properties.  This 
result  is  not  surprising  since  Wallace^5)  found  that  Rec  from  a  quadratic 
tempetaturq  distribution  as  measured  in  his  tests  was  four  times  larger 
than  Rec,  from  an  assumed  linear  Crocco  relation. 


The  rps  deviation  of  the  velocity  data  from  the  theoretical  curve  was  0.8Z, 
which  implies  a^'good  fit".  However,  the  generalized  velocities  for  Y+  <  50, 


* 


which  were  excluded  from  data  fitting  process',  are  considerably  larger  than 
the  theory  reflecting  the  deviation  from  a  linear  profile  that  had  been  ob¬ 
served  near  the  wall.  While  these  data  points  could  be  collapsed  toward 
the  theoretical  curve  by  increasing  uT,  this  would  lead  to  a  large  increase 
in  Cj. 

In  an  effort  to  improve  the  correlation,  the  Van  Driest  formulation  was  modi 
fied  by  replacing  the  Crocco  relation  with  the  experimental  observed  quad¬ 
ratic  law.  (the  assumption  of  constant  pressure  in  Che  boundary  layer  was 
retained) . (21)  The  resulting  curve  fit  was  very  similar  to  that  found  pre¬ 
viously  while  the  new  values  characteristic  parameters  were;  tt  •  1.6,  Cf  m 
4.2  x  10  ,  and  Reg  •  27,000.  It  appears  therefore,  that  the  anomalous 

value  of  rr  and  the  low  value  of  the  fictitious  Re^  can  be  attributed  pri¬ 
marily  to  the  strong  pressure  gradient  normal  to  the  wall. 

A  check  of  the  consistency  of  the  correlation  determined  with  the  assumed 
Crocco  relation  is  shown  in  Figure  32.  While  the  wake  function,  W,  was 
assumed  to  be  given  by  Eqn.  (13),  W  can  also  be  calculated  from  Eqn.  (12) 
after  a  fit  of  the  data  has  been  achieved.  Figure  32,  which  shows  W  de¬ 
termined  from  both  Eqns .  (12)  and  (13)  plotted  against  Y/i ,  indicates  good 
agreement  and  justifies  the  assumption  of  Eqn.  (12). 

v  *  * 

Figure  33  shows  finally  a  plot  of  the  velocity  defect  (ue  -  u  )/u~  versus 
distance  from  the  wall  Y/t  and  illustrates  the  effect  of  the  wake  strength 
parameter  tt.  The  correlation  equation 

(u  *  -  u*)/u_  =  -2.5  In  (Y/5 )  +  1.25(2  -  W)  (14) 

e 

1  ^  . 

which  is  essentially  Eqn.  (13)  written  in  defect  form  with  n  set  equal  to 
0.6,  was  shown  by  Maise  and  McDonald'20'  to  represent  the  incompressible 
measured  profiles  presented  by  Clauser.^2^)  The  agreement  between  Eqn.  (14) 
and  the  measured  velocity  profile,  for  which  tt  *  1.4,  is  quite  poor,  al¬ 
though  it  illustrates  the  lack  of  correlation  that  Maise  and  McDonald^0) 
found  for  non-adiabatic^  boundary  layers.. 


FIGURE  32.  CCMPARSION  OF  COMPUTED  WAKE  FUNCTION 
W,  WITH  ASSUMED  VARIATION,  EQN  (13) 
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SECTION  X 

-A 

SUMMARY  AND  CONCLUSIONS 


flow  properties  in^he^perLnic^^  J0nce™ing  the  measurement  of  mean 
ported  here:  hypersonic  turbulent  boundary  layer  experiments  re- 

(1)  The  results  of  detailed  pitot  pressure  total  r 
static  pressure  surveys  carried  out  9 

demonstrated  the  ln  the  JPL/HWT 

of  the  ceiling  boundary  layer”  H«e»e^  nature 

measurements  vhirh  i ,  r»  c”e  static  pressure 
the  stream^eltrect?™”,?  “8ilgIbl*  variation  in 

niflcant  pressure  gradient  norrsTl'  ^“the  w!"l  I  PrCR»"CI’  3  ' ;  " 
magnitude  of  the  nre«i,ro  al  ff  tn  wal1-  Because  of  »  -ic- 

5°I  of  the  free-stream  pressure^the ’nf1'11  W*S  approxictateiy 
the  boundary  layer  was  taken  i  r*  Pressure  variation  through 
the  data.  ferJl  pressure  Lad,  ^  the  redu«i°"  of 

viously  in  hypersonic  turbulent 'bound^*!  °bserved  pre~ 
walls  as  well  as  on  cones  wed  ^  *ayers  on  wind  tunnel 

have  been  shown  to  increase  wiS^i  f*-t  plates  and»  in  fact, 
gradients,  however,  are  freoue  M  ?reasing  Mach  number.  Such 
precise  measurements  or  In  Z  f°r  lack  of 

perties  are  rel^vely i^e^r  ■  ^  th*  0ther  flofc  pro‘ 

Since  the  observed  prL“e Vari^ns. 

flow  effect  which  is  enhanced  r  k  ^  rePTesent  a  turbulent 
that  further  invesUg^L  esse'tLTf'  5Pe,!ds-  “  *•'  «**» 
tation  of  hypersonic  bo-  ndary  layer  meas  V  accurate  lnterpre- 
the  development  of  impr 

(2>  found  "to 'agree  ™ 

The  quadratic  law  has  been  ob^r^d^"^  ^  b°Undar>'  laTer- 
wall  measurements  and  hat  ho  „°  characterize  nozzle 

to. reflect  the  consequences  of  °th"  -rs 

flPt  plate  data  genera!”  agree  with  t^'cr'’3"51."-  SIn" 
out  the  boundary  layer  it  has  bPPn  H  °  °  theory  through- 

measured  on  nozzle  walls  will  relax  to^he^  ‘f*  profiles 

distances  far  downstream  from  the  throat  1 1**™  relation  at 

that  the  boundary  layer  in  the  JPL/W  \  appears»  therefore, 

librated  in  the  sense  that  rh!  £,7  ^  "0t  COBPle£ely  equi- 

variation  deviates  from  what  is  con^d*1  t®Eperature*velocity 
Plate  flows.  considered  typical  of  flat 

(3)  Using  the  Van  Driest  eenpraii«d  . 

was  correlated  to  the  incom^ressibL^^r  ^h  Wl?city  data 

rms  error  of  less  than  1Z.  PWhile  this  imnf-  he  *ake  *lth  a 
measurements  aeree  wirh  «-k  implies  that  the 

agree  urth  the  general  form  of  the  incompressible 
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profile,  Che  strength  of  the  wake  component  was  twice  as  great 
as  that  usually  found  for  flat  plate  boundary  layers.  This  has 
been  attributed  primarily  to  the  effects  of  the  pressure  gradient 
normal  to  the  wall  since  modification  of  the  correlation  pro¬ 
cedure  to  account  for  the  observed  quadratic  temperature  dis¬ 
tribution  produced  essentially  the  same  results.  In  spite  of 
this  anomaly  the  experimental  skin  friction  coefficient  agreed 
with  the  theoretical  value  within  several  percent. 

(4)  Although  the  mean  flow  measurements  extended  to  0.010  inches 
from  the  wall,  no  indication  of  a  linear  variation  of  velocity 
in  the  sub-layer  was  found.  Possible  causes,  including  the 
effects  of  heat  transfer,  turbulence  and  instrument  errors, 
were  examined  without  conclusive  results.  The  measurements 
indicate  that  the  sub-layer  is  only  0.05  to  0.10  inches  thick, 
which  represents  1  to  2Z  of  the  total  boundary  layer  thickness. 
The  relatively  thin  sub-layer  has  been  attributed  primarily  to 
rhe  high  Reynolds  number.  Re-,  of  the  present  study.  The  sonic 
line  was  found  to  occur  at  Y  *  0.05  inches  so  that  up  to  50Z 
of  the  outer  portion  of  the  sub-layer  is  supersonic. 

(5)  The  stagnation  temperature  profile  calculated  from  hot-wire 
anemometer  data  is  in  good  agreement  with  that  deduced  from 
the  mean  flow  thermocouple  measurements.  A  comparison  of  the 
unit  stream  Reynolds  number  obtained  from  the  hot-wire  and 
mean  flow  measurements  showed  the  same  trend,  but  differences 
as  large  as  507.  were  observed.  The  latter  result  is  not  unex¬ 
pected,  however,  since  it  is  well  known  that  the  hot-wire  pro¬ 
vides  a  relatively  insensitive  measure  of  Reynolds  number. 


o 
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APPENDIX  A 


>  CALCULATION  OF  MEAN  FLOW  PROPERTIES  FROM  HOT-WIRE  DATA 


The  procedure  for  calculating  properties  of  the  mean  flow  in  the  boundary 
layer  from  the  hot-wire  data  was  programed  in  Basic  Language  for  the 
Honeywell  Model  615  Computer  System  in  time-sharing  mode  of  operation  as 
shown  in  the  flow  chart  of  Figure  34.  For  eight  values  of  hot-wire  cur¬ 
rent,  i,  thei corresponding  resistance  was  determined  from  the  measured 
mean  voltage  across  the  wire.  The  method  of  least  squares  was  used  then 
to  provide  a  linear  fit  of  R  versus  i^  from  which  the  slope  :R/of 7  and  Taw 
could  be  found.  Using  the  latter  to  initially  evaluate  the  thermal  con¬ 
ductivity  of  the  fluid,  one  can  obtain  the  zeroth  approximation  to  Nu0 
from  Eqn.  (6)-.  From  NuQ  and  the  calibration  curve  ReQ  versus  NuQ ,  the 
zeroth  approximation  to  ReQ  is  obtained.  The  calibration  curve  of  ", 
versus  ReQ  is  used  then  to  determine  the  recovery  factor  which,  together 
with  T  ,  yields  the  zeroth  approximation  to  the  stagnation  temperature  T,. 
This,  in  turn,  supplies  an  improved  value  of  thermal  conductivity  and  the 
procedure  is  repeated  until  the  successive  approximations  to  the  stagnation 
temperature  converge  to  the  desired  accuracy. 


x 


/ 
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FIGURE  V*  .  CT.OW  CHART  FOR  HOT  WIRF  COMPUTER  PROGRAM 


REFERENCES 


\ 

1.  Kistler,  A.I..,  "Fluctuation  Measurements  in  a  Supersonic  Turbulent 
Boundary  Layer."  Physics  of  Fluids  2.  No.  3,  pp  290-296,  May-June, 
1958. 

2.  Anon.,  "Jet  Propulsion  Laboratory  Wind  Tunnel  Facilities,"  JPL  Tech. 
Memo  33-335,  JPL/CIT  Pasadena,  Calif.,  April  1,  1967. 

3.  J.  Laufer,  H.  Ashkenas  and  A.  Gupta,  private  conmunication. 

4.  Kane,  E.  D. ,  and  Mastlach,  G.  I.,  "Impact-Pressure  Interpretation  in 
a  Rarefied  Gas  at  Supersonic  Speeds,"  NACA  TN  2210,  October,  1950. 

5.  Mathews,  M.  L. ,  "An  Experimental  Investigation  of  Viscous  Effects  on 
Static  and  impact  Pressure  Probes  in  Hypersonic  Flow,"  Hypersonic 
Research^ Memo  44,  C.I.T.,  June,  1958. 

\  1  ' 

6.  Rogers,  K'.  W.,  Wainwright,  J.B.,  and  Touryan,  K.J  ,  "Impact  and  Static 
Pressure  Measurements  in  High  Speed >Flows  with  Transitional  Knudsen 
Numbers,"  pp  151-174,  Supplement  3,  Vol.  2,  Rarefied  Gas< Dynamics . 

edited  by  J.  H.  deLeeuw,  Academic  Press,  New  York,  1966. 

9 

* 

7.  Behrens,  K.,  'Viscous  Interaction  Effects  on  a  Static  Pressure  Probe 
at  M  =  6y"  A1AA  Journal  lt -No.  12,  pp  2864-2866,  December  1963. 

8.  Wagner,  R.D. ,  and  Watson,  R.,  '^Reynolds  Number  Effects  on  the  Induced 

Pressures  of  Cylindrical  Bodies  with  Different  Nose  Shapes  and  Nose 
Drag  Coefficients  in  Helium  at  a  Mach  Number  of  24,"  NASA  TR-R-182, 
Nov.  1963.  .  <* 

,  s 

9.  Fiore,  A.  W.,  "Turbulent  Boundary  Layer  Measurements  at  Hypersonic 
Mach  Numbers,"  ARL  70-0166,  Aerospace  Research  Labs,  l>.  S.  Air  Force, 
Aug.  1970. 

10.  Fisher,  M.  C.,  Maddalon,  D.  V.,  Weinstein,  L.  M.,  and  Wagner,  R.D., 
"Boundary- Layer  Pitot  and  Hot-Wire  Surveys  at  -  20."  AIAA  Journal  9 
No.  5,  pp  826-854,  May  1971. 

11.  Bertram,  M.  H.,  and  Neal,  L. ,  "Recent  Experiments  in  Hypersonic  Tur¬ 
bulent  Boundary  Layers,"  NASA  TMX- 56335,  May  1965. 

/ 

12.  Bushnell,  D.  M.,  Johnson,  C.  B.,  Harvey,  W.  D. ,  ahd  Feller,  W.  V., 
'^Comparison  of  Prediction  Methods  and  Studies  of  Relaxation  in  Hyper¬ 
sonic  Turbulent  Nozzle-Wall  Boundary  Layers,"  NASA  TND-5433,  1969. 

13.  Lee,  R.  E.,  Yanta,  W.  J.,  and  Zeonas,  A.  C.,  'Velocity  Profile,  Skin- 
Friction  Balance  and  Heat  Transfer  Measurements  of  the  Turbulent 
Boundary  Layers  at  Mach  5  and  Zero- Pressure  Gradient,"  NOL  TR  69-106, 
June  16,  1969. 


-68- 


f  REFERENCES  -  (Continued) 

14.  Coles  D.,  "The  Law-of-the-Wak*  in  the  Turbulent  Boundary  Layer  " 
Journal  ot  Fluid  Mechanics  1.  Part  2,  pp  191^226,  1956. 

15.  Walz,  A.,  '^Compressible  Turbulent  Boundary  Layers,"  pp  $00-350 
^*ique^Paris,  j 9^rbule^;  Centre  National  de  Ja  Recherche  Scienti- 

16’  J0bb’  f‘TKV  "inkler’  E*  M*.  and  Persh,  J.,  "Experimental  Investiga'- 

tion  of  Turbulent  Boundary  Layers  in  Hypersonic  Flow,"  NAVORD  REP“3880, 

\ 

17.  Winkler,  E.M.,  and  Cha,  M.  H.,  "Investigation  of  Flat  Plate  Hyper¬ 
sonic  Boundary  Layers  with  Heat  Transfer  at  a  Mach  Number  of  5*2  " 
NAVORD  REP  6631,  1959.  ‘  f' 


18.  Tetervin,  N.,  '<An  Analytical  Investigation  of  the  Flat  Plate^Turbu- 
lent  Boundary  Layer  in  Compressible  Flow,"  NOL  TR  67-39,  May  1967. 

^ '  *e™P’  W’  and  Owens,  F.  K.,  "Nozzle  Wall  Boundary  Layers  at  Mach 

Numbers  20  to  47,”  AIAA  paper  71-161,  1971. 

20.  Maise,  C.,  and  McDonald,  H.,  "Mixing  Length  and  Kinematic  Eddy  Vis¬ 
cosity  ?n  a  Compressible  Boundary  Layer,"  AIAA  Journal  6  No  1 

pp  73-79,  Jan.  1968.  - - *  *  ’ 

21.  R.  Gran,  private  communication. 

22 *  ^Ie*’  D;»  ,,Jhe  Y°ung  Persons  Guide  to  the  Data,"  Proceedings  AFOSR- 

1  Conference  on  Computation  of  Turbulent  Boundary  Layers 

-  1968,  Vol.  II,  compiled  data,  edited  by  D.  E.  Coles  and  E.  A.  Hirst 
Stanford  University,  1969.  ’ 

2j.  Van  Driest,  E.  R.,  "Turbulent  Boundary  Layer  in  Compressible  Fluids," 
Journal  of  Aeronautical  Sciences  18.  No.  3,  pp  145-160,  March  1951. 

24.  Hopkins,  E.  J.,  and  Inouye,  M.,  "An  Evaluation  of  Theories  for  Pre¬ 
dicting  Turbulent  Skin  Friction  and  Heat  transfer  on  Flat  Plates  at 
Supersonic  and  Hypersonic  Mach  Numbers,"  AIAA  Journal  9  No  6 

PP  993-1003,  June  1971.  - - - ’ 

25.  Wallace,  J  E.,  "Hypersonic  Turbulent  Boundary  Layer  Measurements 
Using  an  Electron  Beam,"  AIAAsJournal  1.  No.  4.  nn  7S7-7<iQ  i„r,i 

1969.  also.  IK  CAL  AN-21U-Y-1,  Cornell  Aeronautical  Lb,  Lgu«  1968. 

26.  Clauser ,  F.  H.,  '"The  Turbulent  Boundary  Layer,"  Advances  in  Applied 
Mechanics,  Vol.  1,  Academic  Press,  New  York,  1956~ 


(Page  70  intentiona 1 ly  left  blink) 


4 


•* 

■WKKPHP  FIQl  BUIK-MOT  FIIMED-. 


V 


HMCLASS1FIED 


Sacwoty 


OMkifiMio* 


POPUMCMT  CONTROL  DATA  •  RAD 

W  atilMcl  «<  mmM  mm*'**  m*§ 


i  MtiW  «M  annll  wftft*  «  cl— 


MlfilMTINSACtlVirv  fCaft 

Philco-Ford,  Aetonutronic  Division 
Ford  Road,  Newport  Beach,  California  9^663 


a*  hipokt  ncu*>n  c 

UNCLASSIFIED  ) 


a*  t*Our 


NceotiT  Tin.*  ‘  » 

Mean  Flow  Measurements  in  a  Hypersonic  Solidary  Layer 


«  oescNieTive  »OTi»'cr 


Technical  Report 


_ _ A  \  :  n 

yp»  •*  rmpmrt  an*  mckfira,  *••••>  S. 

ct  October  1970  to  Jul^l" 


iK1971 


I  tUTHOftK  fUM  n— »■  ft«Wft*J 

Laderman,  Arnold  J.  and  Demetriades,  Anthony 


MPONT  DATE 

August  1971 


a  CONTRACT  o*  *»»**  MO 

F04701-70-C-0130 

ft  MSJIC'NO  _ 


an  TOTAL  NO ,sr 
82 


76  MO  or  RIM 
26 


•  a  OKI  4  UR  A  TO  MtfOMT  NUMOIVij 

TBD  *  t 


othi*  »iro«T 
ftii  now 

U-4950 


»  M 


io  *v*ii-*Bu.iTv/u««iT*TiooMioncM  Each  transmittal  of  this  document  outside  the  Depart¬ 
ment  of  Defense  most  have  prior  approval  of^AMSO  (RNSE)  Norton  Air  Force  Base 
California  92409.  The  distribution  of  this  report  is  limited  because  it  contai 
technology  requiring  disclosure  only  within  the  Department  of  Defense 


1 1  SO0M»l  VMCMTftftV  NOTH 


ia  •■►OMOIUII8  MiuTftftv  activity 


Advanced  Research  Projects  Agency 


*  ^"experimental  investigation  was  carried  out  to  determine/the  structure  of  the 
turbulent  boundary  layer  on  the  wall  of  the  Jet  Propul sW  Laboratory  Hypersonic 
Wihd  Tunnel  at  a  free-str earn  Mach  number  rf  9.3/.  Profiles  o  . 

were  obtained  from  pitot  pressure,  static  pressure,  and  total  temperature  survey, 
made  through  the  6  inch  thick  boundary  layer  at  a  station  a60  inches  from  the 
nozzle  throat.  Tests  were  conducted  primarily  for  the  following  con  *  . 

free-stream  Reynolds  number  of.  127, OCX)  per  inch,  corresponding  to  a  Reynolds 
number  based  on  momentum  thickness  of  36,800,  and  ^  wall  to  free-stream  total 
temperature  ratio  of  0.383.  A  cursory  examination  of  the  boundary  layer  wa8  aV  j 
made  at  a  unit  stream  Reynolds  number  of  67,000.  The  stat;i1P^ssur*/",^“r^T 
although  indicating  negligible  streanwise  variation,  revealed  the 
significant  pressure  gradient  normal  to  the  wall,  with  the  pressure  at  the  “a 
approximately  457,  greater  than  its  free-stream  value.  The  data  1 

profile  ot  total  temperature  ratio  vefsus  velocity  ratio  agrees  with  the  linfear 
Crocco  theory  in  thesub- layer  region,  but  in  the  outer  portion  of  ^boundary 
layer,  closely  follows  the  quadratic  relation  which  has  been  observed  to  charac 
tSU  »zzUy«ll  neasurenents .  The  Indicate  that  the  »ub-l.,yer*s 

0.05  to  0.10  inches  thick  which  corresponds  to  only  seyeral  percent  of  the  total 
boundary  layer  thickness.  Correlation  of  the  velocity  profile  with  the  conven¬ 
tional  incompressible  profile  shows  poor  agreement  in  the  sub-layer  and  wall-of 
the -wake  regions ,  with  the  latter  attributed  primarily  to  the  pressure  gradient 
normal  to  the  wall.  However,  the  experimental  skin  friction  coefficient  g 
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with  the  value  predicted  by  the  Van  Driest  theory  within  several  percent. 
Finally,  the  total  temperature  profile  determined  from  a  hot-wire  anemometer 
traverse  througi  the  boundary  layer  is  in  good  agreement  with  the  thermocouple 
measurements . 


